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ABSTRACT 

Over-the-horizon  radar  (OTHR)  and  microwave  radar  networks  can  to¬ 
gether  generate  track  data  over  a  wide  surveillance  region.  However  the  data 
is  often  subject  to  ambiguity  and  uncertainty  due  to  the  complexities  of  the  HF 
signal  propagation  environment,  which  give  rise  to  multipath  OTHR  tracks,  as 
well  as  ambiguities  in  target  associations  between  multiple  microwave  radars. 
This  report  describes  an  association  and  fusion  algorithm  which  deals  with 
both  sources  of  uncertainty.  The  algorithm  is  capable  of  fusing  OTHR  multi- 
path  tracks  and  non-OTHR  tracks  (e.g.  microwave  radar  or  GPS),  as  well  as 
dealing  with  multipath  tracks  from  OTHR  networks.  The  algorithm  achieves 
this  through  a  very  general,  model  based,  approach  which  deals  with  multi- 
path  effects  as  well  as  asjmchronicity  between  sources  of  data.  Importantly, 
the  approach  incorporates  track  history  in  its  computation  of  association  prob¬ 
abilities  and  fused  estimate  calculations,  thus  exploiting  temporal  as  well  as 
instantaneous  spatial  relationships  between  tracks. 
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Over-the-Horizon  Radar  Multipath  and  Multisensor  Track 
Fusion  Algorithm  Development 


EXECUTIVE  SUMMARY 

The  use  of  over-the-horizon  radar  (OTHR)  in  conjunction  with  microwave  radar  (MR) 
networks  enables  surveillance  over  a  large  region  of  land /sea,  with  potentially  greater  re¬ 
liability  and  accuracy  in  target  acquisition  and  tracking  over  a  single  radar.  For  example, 
the  accuracy  and  clarity  of  OTHR  data  can  be  enhanced  by  fusing  OTHR  tracks  with  asso¬ 
ciated  microwave  radar  tracks  which  are  known  with  greater  precision  and  axe  not  subject 
to  multipath  effects  of  the  sort  that  plague  OTHR.  To  achieve  this  potential,  substantial 
challenges  which  occur  due  to  ambiguities  and  uncertainties  in  target  associations  between 
multiple  microwave  radars  as  well  as  multipath  OTHR  tracks,  must  be  overcome.  Multiple 
target  tracks  arise  in  OTHR  surveillance  due  to  multiple  ionospheric  propagation  paths 
between  targets  and  radar  locations.  The  use  of  multiple  microwave  radars  also  leads  to 
multiple  tracks  for  single  targets.  Hence  target  track  association  and  fusion  algorithms  for 
OTHR/MR  networks  are  required  to  give  a  unified  picture  of  the  surveillance  region  with 
one  track  per  target  in  ground  coordinates. 

This  report  describes  an  algorithm  for  the  association  and  fusion  of  multipath  OTHR 
tracks  as  well  as  microwave  radar  tracks  which  fall  within  the  OTHR  surveillance  region. 
Both  multiple  microwave  radars  and  multiple  OTHRs  can  be  accommodated  by  the  algo¬ 
rithm.  Global  Positioning  System  (GPS)  reports  from  commercial  aircraft,  air  lanes  and 
airfield  locations  could  also  be  incorporated  if  required.  This  algorithm  can  be  used  to 
enhance  the  value  of  the  OTHR  surveillance  data  from  the  Jindalee  Facility  Alice  Springs 
(JFAS)  and  the  Jindalee  Over-the-horizon  Radar  Network  (JORN)  for  air  picture  compi¬ 
lation  and  to  aid  the  integration  of  OTHR  into  a  multi-sensor  air  surveillance  system.  At 
the  time  of  writing  this  report,  testing  of  the  algorithm  with  multipath  OTHR  tracks  has 
commenced  on  a  test-bed  as  well  as  on  a  prototype  version  in  an  operational  OTHR.  The 
results  of  testing  to  date  indicate  that,  at  the  present  state  of  development,  the  algorithm 
can  be  used  in  a  semi-automated  fashion  to  give  advice  to  an  operator.  However,  higher 
levels  of  automation  can  be  expected  with  further  development. 
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all  the  (radar  coordinate)  track  data  available  up  to  and  including  update  k 
£'{.}  mathematical  expectation 
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(fc  —  1)  and  k 
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matrix  P 
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Pj"-’  {k)  covariance  matrix  of  yj*-’  (fc),  ie,  P{"^  (fc)  =  E  {k)  y”*^  (A:)^| 
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(fc)  the  estimates  xjji  (k) ,  (fc) 
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=  ^;{C'w(*rw)'} 


Sj  the  set  {rrij  :  mj  —  rrij-i  if  rij  =  nj-i,  ruj  =  mj-2  if  nj  =  nj-2,.-.,  ruj  =  mi  if 
rij  =  ni} 

(k)  covariance  matrix  of  (k)  —  z^  {k)J 
ti  target  i,  i  =  1,  ..,r 

T  the  number  of  targets  associated  with  a  hypothesis 

rJ”-’  (fe)  covariance  matrix  of  (A:)  —  (A)  j 

Tj  track  j,  j  —  in  radar  coordinates 

the  event  “track  Tj,  j  —  1, ..,  J  is  associated  with  propagation  via  propagation  path 
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(A)  cross  covariance  matrix  of  (A)  and  z^’  (A),  ie,  (A)  =  E  (A)  ^zj'^  (A)^  ^ 

(A)  cross  covariance  matrix  of  z^^  (A)  and  xf  (A),  ie,  (A)  =  E  (A)  (A))^| 

V  (A)  zero-mean  white  Gaussian  target  process  noise  at  update  (A) 

Vm  “volume”  in  measurement  space  in  which  a  target  can  be 

Vs  “volume”  in  state  space  in  which  a  target  can  be 

Xi  (A)  the  state  vector  for  target  ti  (ie,  its  true  state)  at  time  A 
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state  estimates  at  time  A  and  the  fused  estimate  from  time  A  —  1 

x^^  (A)  error  corresponding  to  the  estimate  x^  (A),  ie,  x^  (A)  =  Xj  (A)  —  x^  (A) 

(A)  the  estimate  in  ground  coordinates  corresponding  to  xpj  (A)  assuming  propagation 
path  rrij 
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(k)  error  corresponding  to  the  estimate  (k),  ie,  y^  ■'  (fc)  —  Xj  (k)  —  yj  ^  (k) 

zj’’^  (/c)  the  measurement  in  ground  coordinates  corresponding  to  Sj  (k)  assuming  propa¬ 
gation  path  rrij 

(fc)  the  measurement  prediction  (the  expected  value)  for  target  ti  at  time  k  resulting 
from  the  fusion  of  h  earlier  state  estimates  at  time  k  and  the  fused  estimate  from 
time  A:  —  1 

(k)  error  corresponding  to  the  measurement  Zj  ^  {k),  ie,  z^  ^  (k)  =  Hxi  (k)  —  Zj  ■’  (A:) 


DSTO-RR-0223 


1  Introduction 

In  skywave  over-the-horizon  radar  (OTHR),  ionospheric  refraction  of  HF  signals  is  ex¬ 
ploited  to  detect  and  track  targets  well  beyond  the  Une-of-sight  horizon,  thus  providing 
surveillance  over  much  larger  areas  than  is  possible  with  conventional  ground  based  mi¬ 
crowave  radar.  Unfortunately,  ionospheric  propagation  conditions  are  usually  such  that 
several  propagation  paths  exist  between  the  target  and  the  radar  site,  thus  giving  multiple 
resolved  measurements  for  a  single  target.  Figure  1  shows  a  simple  model  of  the  multiple 
propagation  path  phenomenon.  The  figure  shows  a  transmitter  at  point  O  and  target  at 
point  P,  with  two  ionospheric  layers,  the  E  layer  and  the  F  layer,  via  which  the  transmitted 
and  returned  signals  can  travel.  A  signed  transmitted  from  point  O  can  travel  via  one  of 
two  paths  to  point  P,  and  similarly,  a  signal  returning  from  point  P  to  point  O  can  travel 
by  one  of  two  paths.  If  the  path  lengths  for  the  four  path  combinations  are  sufficiently 
different  for  resolved  returns  to  occur,  this  leads  to  four  “apparent”  targets  when  there  is 
actually  only  one. 


Figure  1:  Multipath  propagation  via  ionospheric  E  and  F  layers 

Cmrent  practice  is  for  OTHR  tracking  to  be  implemented  in  radar  coordinates  [group 
range,  group  range  rate,  apparent  azimuth  and  apparent  azimuth  rate),  with  the  major 
reason  being  to  limit  computational  requirements.  The  result  is  that  when  multipath 
propagation  conditions  exist,  multiple  radar  tracks  are  formed  for  each  target.  In  situ¬ 
ations  where  there  is  an  unknown  number  of  targets,  the  challenge  is  then  to  correctly 
associate  radar  tracks  with  targets  and  to  subsequently  estimate  target  locations  in  ground 
coordinates  by  fusing  associated  radar  track  data.  Many  plausible  track-to-target  associ¬ 
ations  may  exist,  so  the  association  ambiguity  must  be  resolved. 

There  have  been  several  research  efforts  investigating  the  association  and  fusion  of 
multipath  OTHR  tracks  which  have  been  reported  in  the  literature  [11],  [12],  [8],  [28]  as 
well  as  an  approach  which  directly  performs  tracking  of  a  target  in  ground  coordinates 
using  multipath  measurements  in  the  presence  of  clutter  [22].  Of  these,  the  work  de¬ 
scribed  in  [8],  [28],  [22]  represent  early  efforts  sponsored  by,  or  performed  in.  Surveillance 
Systems  Division  (SSD).  Subsequently,  a  substantial  effort  was  made  in  SSD  to  develop  a 


1 


DSTO-RR-0223 


new  algorithm  for  multipath  track  fusion  (MPTF)  by  Percival  and  White.  In  the  MPTF 
approach,  all  track-to-target  association  hypotheses  are  recursively  constructed,  the  prob¬ 
ability  of  each  hypothesis  is  calculated  and  the  fused  target  states  within  each  hypothesis 
are  evaluated.  The  approach  is  described  in  references  [17],  [18]  and  [19]. 

Although  the  MPTF  technique  was  a  promising  approach  to  the  problem,  it  suffered 
from  the  following  deficiencies.  The  algorithm  treats  each  track  update  independently, 
which  results  in  fusion  of  track  estimates  without  taking  into  account  track  history.  Addi- 
tionaly,  the  algorithm  does  not  have  an  effective  hypothesis  management  strategy,  making 
it  impractical  in  all  but  the  most  trivial  circumstances.  In  the  absence  of  a  hypothesis 
management  strategy,  the  number  of  association  hypotheses  produced  by  the  MPTF  al¬ 
gorithm  quickly  becomes  too  large  for  real  time  execution.  Some  preliminary  work  on 
overcoming  the  difficulties  in  hypothesis  management,  which  was  done  by  Percival  and 
White,  is  described  in  Refs.  [18]  and  [19];  however  the  approach  proved  ineffective  for  rea¬ 
sons  that  will  be  described  later.  A  further,  and  very  restrictive,  limitation  of  the  MPTF 
algorithm  is  its  inability  to  deal  with  asynchronous  track  data.  The  outcome  is  that  track 
data  from  multiple  asynchronous  sensors  cannot  be  accomodated.  Thus  the  algorithm 
cannot  exploit  synergies  that  exist  when  multiple  OTHRs,  microwave  radars  and  other 
sensors  are  used  in  conjunction. 

The  use  of  multiple  OTHRs  in  conjunction  with  microwave  radar  networks,  as  well  as 
other  track  level  information  such  as  GPS  reports,  offers  the  potential  of  surveillance  over 
a  large  region  of  the  earth’s  surface,  with  substantially  improved  reliability  and  accuracy  in 
target  acquisition  and  tracking  over  a  single  OTHR.  Skywave  OTHR  offers  the  advantage 
of  radar  coverage  over  a  large  area  of  the  earth’s  surface  as  well  as  the  ability  to  detect 
low  flying  targets  well  past  the  horizon,  but  suffers  from  drawbacks  which  include  multi¬ 
ple  resolved  measurements  of  individual  targets  and  generally  inaccurate  registration.  In 
contrast,  microwave  radar  offers  much  higher  measurement  accuracy  and  is  largely  free  of 
the  multiple  measurement  problem  associated  with  OTHR.  However  a  single  microwave 
radar  has  a  very  limited  coverage  area  in  comparison  with  OTHR,  particularly  in  the 
case  of  a  ground  based  radar  and  low  flying  targets.  GPS  reports  have  the  advantage  of 
even  greater  accuracy,  but  have  only  limited  availability  and  are  certainly  not  available 
for  unfriendly  targets.  Fusing  the  OTHR  measurements  with  microwave  radar  measure¬ 
ments  and  GPS  reports,  where  available,  thus  enables  the  advantages  of  each  sensor  to  be 
utilized  in  areas  of  overlapping  coverage,  to  achieve  better  accuracy  and  reduced  associ¬ 
ation  ambiguity.  Also,  a  target  that  is  temporarily  in  a  Doppler  blind-zone  of  one  radar 
may  be  visible  to  another  radar;  hence  track  continuity  can  be  enhanced  by  multi-sensor 
fusion.  Furthermore,  by  taking  advantage  of  spatial  continuity  of  the  ionospheric  layers, 
information  gained  by  fusing  OTHR,  microwave  radar  and  GPS  measurements  in  areas  of 
overlapping  coverage  can  be  used  to  more  accurately  determine  ionospheric  layer  heights 
in  these  areas,  and  by  extrapolation,  in  surrounding  areas.  Hence  the  performance  of  the 
OTHR  itself  can  be  improved,  not  only  where  there  is  overlapping  coverage,  but  also  for 
arefis  nearby. 

To  perform  fusion  of  OTHR,  microwave  radar  and  GPS  measurements,  the  major  task 
is  again  to  correctly  associate  measurements,  which  in  this  case  are  from  multiple  sensors, 
with  targets  and  to  subsequently  estimate  target  locations  in  ground  coordinates  by  fusing 
associated  radar  (and  GPS)  track  data  to  give  a  unified  picture  of  the  surveillance  region. 
Because  of  the  similarities  between  the  tasks  of  fusing  multipath  tracks  to  that  of  fusing 
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multisensor  tracks,  it  seems  natural,  as  well  as  very  desirable,  that  a  single  algorithm 
should  be  developed  for  both  problems. 

This  report  describes,  a  new  dynamic  multipath  track  fusion  (DMPTF)  algorithm 
which  is  capable  of  performing  association  and  fusion  of  multipath  OTHR  traxdcs  from 
multiple  OTHRs  as  well  as  multiple  microwave  radar  tracks  which  fall  within  the  OTHR 
surveillance  region.  GPS  reports  from  commercial  aircraft  can  also  be  incorporated  if 
required.  The  new  algorithm  builds  on  the  earlier  work  done  in  SSD  by  Percival  and 
White  [17],  [18],  [19]  by  providing  major  extensions  which  overcome  the  limitations  of 
the  earlier  algorithm.  The  DMPTF  algorithm  treats  all  the  sources  of  track  data  in  the 
same  fashion,  leading  to  a  simple  and  elegant  set  of  association  and  fusion  equations. 
This  algorithm  can  be  used  to  enhance  the  value  of  the  OTHR  surveillance  data  from  the 
Jindalee  Facility  Alice  Springs  (JFAS)  and  the  Jindalee  Over-the-horizon  Radar  Network 
(JORN)  for  air  picture  compilation,  to  aid  the  integration  of  OTHR  into  a  multi-sensor 
air  smveillance  system,  and  thereby  contribute  to  the  command  and  control  capabilities 
of  the  Australian  Defence  Force  (ADF).  Early  trials  of  the  algorithm  have  been  performed 
on  a  test-bed  as  well  as  on  an  operational  radar,  giving  very  promising  results. 

The  report  describes  the  advances  made  to  the  earlier  MPTF  approach,  the  key  ad¬ 
vances  being: 

•  an  algorithm  which  enables  temporal  as  well  as  spatial  track  behaviour  to  be  ac¬ 
counted  for  during  the  association  and  fusion  processes, 

•  a  robust  strategy  for  the  pruning  of  association  hypotheses,  which  enables  the  effi¬ 
cient  implementation  of  the  algorithm  for  a  much  larger  number  of  tracks  than  was 
previously  possible, 

•  the  development  of  an  algorithm  for  dealing  with  changes  in  the  number  of  tracks 
during  consecutive  updates  (a  common  phenomenon  with  multipath  tracks), 

•  equations  for  using  measurements  in  place  of  track  estimates  in  the  association  and 
fusion  processes,  the  aim  of  these  equations  being  to  enable  the  reduction  of  effects 
due  to  track  dependencies, 

•  the  development  of  an  algorithm  that  can  be  used  to  associate  and  fuse  asynchronous 
non-OTHR  tracks  (eg.,  microwave  radar  or  GPS)  with  the  multipath  OTHR  tracks 
of  multiple  OTHRs  in  a  single  fusion  process. 

In  the  earlier  work  by  Percival  and  White,  the  association  and  fusion  equations  are 
obtained  by  analogy  to  multisensor  work  described  in  the  literature.  In  addition  to  the 
advances  summarized  in  the  previous  paragraph,  a  further  contribution  of  this  report  is 
the  first  principles  derivation  of  the  original  MPTF  algorithm,  which  has  confirmed  errors 
in  the  equations  derived  using  the  analogy  approach  of  the  earlier  work.  Furthermore,  the 
first  principles  approach  is  subsequently  used  to  provide  a  rigorous  derivation  of  the  new 
dynamic  algorithm  as  well  as  equations  for  using  measurements  in  place  of  track  estimates 
in  the  association  and  fusion  processes.  Another  contribution  of  the  report  is  a  literature 
survey  and  discussion  of  possible  approaches  for  dealing  with  track  dependence. 
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This  report  is  organized  as  follows:  In  Section  2,  a  first  principles  derivation  of  the 
earlier  “static”  MPTF  algorithm  is  given.  Section  3  discusses  issues  regarding  the  com¬ 
putational  complexity  of  the  MPTF  algorithm  and  gives  reasons  why  hypothesis  pruning 
is  necessary  for  the  algorithm  to  be  able  to  be  implemented.  In  section  4,  the  develop¬ 
ment  of  a  hypothesis  pruning  algorithm  is  described,  and  subsequently,  in  section  5,  its 
implementation  is  discussed.  Section  6  presents  the  new  dynamic  multipath  track  fusion 
(DMPTF)  algorithm,  and  investigates  approaches  for  dealing  with  track  dependence.  Sec¬ 
tion  7  then  goes  on  to  describe  the  results  of  initial  performance  testing  and  assessment. 
In  section  8,  a  modified  version  of  DMPTF  is  derived  in  which  measurements  are  used  in 
the  association  and  fusion  process  in  place  of  track  estimates.  Section  9  describes  how  the 
DMPTF  algorithm  can  be  used  for  fusion  of  tracks  from  multiple  sensors  such  as  multiple 
OTHRs  or  OTHRs  and  microwave  radars.  Finally,  in  section  10  a  summary  of  the  work 
and  conclusions  are  presented  as  well  as  an  outline  of  future  work  that  is  recommended. 
The  reader  is  also  referred  to  [24]  and  [25],  which  present  some  of  the  algorithms  and 
results  described  in  this  report,  but  in  substantially  less  detail  than  is  given  here. 
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2  Static  Multipath  Track  Fusion 

The  static  MPTF  algorithm  [17],  [18],  [19]  generates  association  hypotheses  by  recur¬ 
sively  constructing  a  hypothesis  tree  at  each  update  of  the  tracks.  Each  path  through  the 
tree  corresponds  to  a  unique  association  of  radar  tracks  with  targets  and  available  prop¬ 
agation  paths.  The  algorithm  evaluates  the  probability  of  each  hypothesis  and  computes 
estimates  of  the  target  states  within  each  hypothesis.  Key  points  to  note  are  that  the 
tracks  must  be  updated  simultaneously  and  each  update  is  treated  independently  by  the 
MPTF  algorithm.  A  description  of  the  static  algorithm  is  given  in  references  [17],  [18]  and 
[19];  however,  the  earlier  version  was  developed  by  analogy  with  other  work  and  was  not 
rigorously  derived.  There  are  also  some  mathematical  errors  in  the  earlier  equations.  For 
these  reasons,  and  because  the  static  algorithm  forms  the  basis  of  the  initial  creation  of  the 
hypothesis  tree  in  the  DMPTF  algorithm,  it  was  necessary  to  derive  the  static  equations 
rigorously  from  first  principles.  The  derivation  is  presented  in  this  section. 


Figure  2:  Four  multipath  tracks  in  radar  coordinates 

Consider  a  collection  of  tracks  tj,  j  —  1, ...,  J  in  radar  coordinates.  Figure  2  shows  an 
example  with  four  tracks  (ie,  J  =  4),  which  often  occurs  with  a  single  target.  Let  us  first 
consider  update  (time)  k  =  0.  At  time  k  —  0  the  J  tracks  are  represented  by  the  estimates 
V’j(O),  j  =  Given  prior  information  on  the  propagation  path  transformations 

from  radar  to  ground  coordinates  [17],  the  corresponding  estimates  in  ground  coordinates 

(0)  and  their  covariances  (0)  can  be  calculated  for  each  possible  propagation  path 
mj.  Using  this  information,  a  hypothesis  tree  for  the  J  tracks  is  created  recursively, 
starting  from  the  root  and  building  the  tree  by  considering  one  track  estimate  at  a  time 
until  all  the  estimates  have  been  processed.  The  hypothesis  tree  is  used  as  a  tool  to 
visualize  all  the  target  and  path  associations  that  are  possible  for  the  J  tracks.  Note  that 
for  the  remainder  of  this  section  the  time  index,  (0),  is  not  written,  for  the  sake  of  brevity. 

To  begin,  let  Xnin2-.nK  denote  the  composite  hypothesis  that  Tj  is  associated  with 
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target  n^,  j  =  K  <  J,  ie,  ri  is  associated  with  ni,  T2  is  associated  with 

W2r")  '^K  is  associated  with  uk-  This  type  of  hypothesis  will  be  referred  to  as  a  path 
independent  hypothesis  in  the  remainder  of  this  report.  As  an  example,  figure  3  shows 
a  path-independent  hypothesis  tree  for  a  cluster  of  three  tracks  {J  =  3).  The  first  track 
is  assigned  to  target  ti  yielding  hypothesis  Ai.  The  second  track  can  be  due  to  the  first 
target,  ti,  or  a  second  target,  t2,  giving  the  set  of  hypotheses  {An,  A12}.  If  the  first  two 
tracks  were  due  to  target  ti ,  the  third  track  could  be  due  to  ti  or  could  be  due  to  a  new 
target  t2,  giving  the  hypotheses  (Am,  A112}.  If  the  first  two  tracks  are  due  to  two  different 
targets,  that  is  ti  and  t2  respectively,  then  the  third  track  can  be  due  to  either  ti,  t2  or 
a  new  target  ts  giving  the  hypotheses  {A121,  A122,  Ai23}.  Thus  for  three  tracks  there  are 
five  possible  path  independent  hypotheses  (Am,  Am,  A121,  A122,  A123}.  Note  that  in  the 
general  case,  assuming  that  there  are  no  false  tracks,  the  values  that  rij  can  take  lie  in 
the  set  {1, 2, ...,  }>  where  is  the  number  of  path  independent  hypotheses 

that  originate  from  hypothesis  For  example,  consider  A12  in  figure  3,  where 

Bi2  -  3. 


Figure  3:  Path  independent  hypothesis  tree  for  three  tracks. 

Now  let  denote  the  composite  hypothesis  that  tj  is  associated  with  target  nj 

and  propagation  path  mj,  j  —  1,2,...,K,  K  <  J.  This  type  of  hypothesis  wOl  be  referred 
to  as  a  path  dependent  hypothesis.  The  values  that  can  be  taken  by  nj  are  as  before,  and  ruj 
lies  in  the  set  (1,2, ...,  Mj}  where  Mj  is  the  number  of  propagation  paths  associated  with 
'tpj.  Note  that  Mj  is  known  a  priori,  as  it  is  supplied  as  part  of  the  coordinate  registration 
(CR)  information.  If  it  is  assumed  that  no  two  resolved  ground  tracks  that  are  due  to 
the  same  target  can  be  associated  with  the  same  propagation  path,  then  the  values  that 
nij  can  take  lie  in  the  set  (1, 2, ...,  A/j}  n  Sj  where  Sj  ={mj  ;  mj  =  mj_i  if  Uj  —  Uj-i, 
mj  —  mj-2  if  nj  =  nj_2,...,  mj  =  mi  if  Uj  =  ni}  and  Sj  is  the  complementary  set  of  Sj. 
This  assumption  is  made  for  the  remainder  of  the  report.  Clearly,  each  path-independent 
hypothesis  Xni,..,nj  corresponds  to  a  set  of  path- dependent  hypotheses  }.  The 

relationship  between  path-dependent  and  path-independent  hypotheses  can  be  illustrated 
using  the  previous  example  of  a  cluster  of  J  =  3  radar  tracks  and  letting  the  number  of 
propagation  paths  be  Mj  =  2.  Consider  the  sequence  of  hypotheses,  Ai  — +  A12  — ♦  A121,  in 
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the  path  independent  hypothesis  tree  of  Figure  3.  This  sequence  can  be  decomposed  into 
the  path  dependent  hypotheses  shown  in  Figure  4. 


Figure  4:  Path  dependent  hypotheses  corresponding  to  path  Ai  — >  A12  — >  A121  of  the  path 
independent  hypothesis  tree. 

Consider  now,  in  detail,  the  creation  of  the  hypothesis  tree  and  the  calculation  of  the 
hypothesis  probabilities  and  fused  estimates. 

Firstly,  let  us  define  the  following: 

=  the  event  “track  Tj,  j  =  is  associated  with  propagation  via  propagation 

path  mj,  nrij  e  {1,2, M^}”. 

0”^  =  the  event  that  “track  Tj,  j  =  1, J  is  associated  with  target  nj,  nj  €  (1, 2, J}”. 

From  the  first  estimate,  ipi,  the  propagation  path  transformations  can  be  used  to 
calculate  and  its  covariance  for  each  propagation  path  mi  =  1,  Under  the 

assumption  that  there  are  no  false  tracks,  there  is  only  one  target  association  possible  for 
rpi,  namely  target  1.  The  probabihty  of  each  of  the  path  dependent  hypotheses  is  then 

P{>^r\i’i}  =  P{eT\4>\M  mi  =  l,..,Mi 

where 

P  =  1 

hence 

=  p{^rkUi} 

=  p{0f'} 

Let  be  the  prior  probability  of  ;  this  is  estimated  using  physical  measurements  of 
the  ionosphere.  Hence  P  =  /3^‘,  and 

P{XT'\^i}=^f3r  mi  =  l,..,Mi  (1) 
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For  each  subsequent  estimate,  ijjj,  j  =  the  recursive  equation  for  calculating 

the  probability  of  from  be  derived  as  follows. 

Consider  the  second  estimate,  V’2;  using  the  propagation  path  transformations  we  can 
then  calculate  and  its  covariance  for  each  propagation  path,  ie,  m2  =  1,..,M2, 
where  M2  is  the  number  possible  paths  for  T2-  Define  now  =  the  (single  path  depen¬ 

dent)  hypothesis  that  Tj  is  associated  with  target  rij  €  {l,  2, }  and  propagation 
path  rrij  G  {1, 2, Mj}. 

Consider  now  hypothesis  track  r2  is  associated  with  target  712  and  with 

propagation  path  m2.  The  probability  of  the  hypothesis  is 


p  \Mi  ,xr}  =  P  {er .  1^2,  ,  A^  ^ } 

2  M2 
n2=l  m2  =  l 

p {<>2  AT)P {OT  (■»?  \i>uK'} 

2  Afo 
n2=l  m2=l 


Now,  in  the  above  P  {0^'^  ,  A^‘ }  =  P  {0'^^}  =  0^^ ,  m2  =  1,..,M2  because 

the  conditions  give  no  information  regarding  the  probability  of  6^^.  With  regard  to 
p  {4>T  IV'I)  '^■T^  }  >  ^2  =  1)2,  there  are  (at  least)  two  approaches  that  may  be  considered: 

Approach  1:  Assume  that  we  know  nothing  about  the  target  density,  and  simply 
assume  that  both  4>\  and  4>2  are  equally  likely,  ie,  P  ]  =  ;^  =  ^  and 

P  {(^2  IV’l)  A^^'  }  ~  “  2  •  "^bis  is  the  approach  that  will  be  developed  here. 

Approach  2:  Assume  that  we  know  something  about  the  target  density  a  priori  and 
use  this  to  determine  P  {(p2^  \ipiAT^  }■  This  approach  will  not  be  investigated  in  this 
report. 

Now,  in  the  implementation  of  MPTF  the  ground  co-ordinate  estimates  yj*-’  are  used 
to  perform  the  probability  calculations,  hence  we  need  an  equation  in  terms  of  them.  To 
obtain  this,  wherever  we  encounter  the  conditioning  event  we  can  replace  Tpj  with 

,  ie,  replace  the  estimate  in  radar  coordinates  with  the  equivalent  estimate  in  ground 
coordinates  given  the  stated  propagation  path.  Hence  we  can  easily  obtain 


pfA™iv>2,*.Ar}  = 

p (vT  ■K')P (g’  !<»”■  !<r‘  AT)P{<l’V\yT'. K' ) 

2  Mo 

E  E  p(yri&r-0^^yr,K')P{^rhr-yr^^T'}PW,lyr^K'} 

n2—l  m2—l 
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Substituting 


into  the  above  gives 


P{^K^\MuK'}  = 


p{yTK\'H\yT,K')l^ 


m2 


E"’=.  p  (yP  \«P  ,4\vT\K')iV^ 

The  probability  of  the  composite  hypothesis  is  hence  given  by 

piK^Mi}  =  pfA-/,AriV'2,V’i} 


pjyr  l^^.<^2^y^^A^)/3^ 


m2 


Ei=i  Ezu  p  {yp  \or ,  >  yr .  Ar  ^ 


m\ 


where  it  has  been  assumed  that 

p{AriV’2,V'i}=i’{AriV'i} 

ie,  that  the  probability  of  is  independent  of  '\p2  prior  to  any  hypothesis  being  made 
about  the  target  and  propagation  path  associated  with  ^’2- 

Consider  now  the  general  case  of  calculating  the  probability  of  from  An/.;;I^^7* . 

We  then  have  the  following 

nij  €  {1, 2, Mj}  ,  Tij  €  'll,  2, J  ~  Ij  ■•)  ^ 

Let  us  now  introduce  the  following  shorthand  notation:  =  ^1, xpK  >  where  again  the 

time  index,  (0) ,  is  omitted  for  brevity.  Then 

p (fe ) p  W'.’iV  !*'"-■, ASVnT:' } 


p  {<Pk  IST .  •i’W  AS' :„Ti‘ )  X 

{<"  }  X  p  {r/ 1*''-'. 


^mi  ..iriK-i 


■■rriK-l  \ 

riK-l  J  J 


Bni..TiK^_l  Mk 

E  E  p  (fe  \0T ,  .  aS‘.;„Tt'  )p{epA7\^'^-',  C.inT-i' } 

nff=l  771^  =  1 

p{9r  }  X  F{i.7  )  J 


E  E 

71^=1  mK=l 


p(fe|«5x,4*.'l'''"'.AM-T‘)x 

[  P  {P™  l«7,i''f-‘,AS‘:„T-T‘ }  P  {C.  |»'''"‘,aS':«Tt'  } 
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As  stated  earlier,  the  relevant  estimates  for  calculating  the  probabilities  are  the  ground 
co-ordinate  estimates  .  The  equation  in  terms  of  them  is 


^ni..nK_i 


} 


\p{yr  i»r.( 
xP{erW 

Pk  >2/1  '"'yK-l  ''^ni-.TiK-i  J 

^  iPl  '"'yK-l  ''^ni-.nK-i  } 

1  mi  %mi..mK-_i  1 

2/1  '"'yK-1  )Ani..nK-i  / 

Mk 

E  E 

nK=l  rnK  =  l 

PyyK  rK  'Vk  >yi  '■•lyK-l  i^ni..nK-\  / 

^P  K  \Pk  >2/1  '"'yK-1  '  S 

D  f  1.,  rui  ruK'-l  xmi-.m/c-i  1 

xP  ,  |yj  ^-lyK-l  iAni..nK-i  /  J 

Substituting 


p{e^^ 


4>' 


PK 
K  > 


ruK-i 

’VK-I 


1 


} 


m/c-i  xTUi.-niK-i  \ 
Vk-X  J 


/3 


rriK 

K 


1 


B 


m.-TiK-i 


into  the  above  gives 


« mi-.m/f-i 
Ani..nK--i 


„  >TiK-i  \»ni..mK-i  \  ^rriK 

yK-1  )  Pk 


^"1  ■•"K'-l  Mk  _  .  . 

E  E  .i/r 

nK  =  l  ”IK-  =  1 


niK-i 

yK-1 


)/?: 


<mK 

K 


The  probability  of  the  composite  hypothesis  P  }  is 


p  ,  A;r;.;nT-T‘  } 

Now,  applying  conditional  probability  and  assuming  that^ 


■r  •[Ani..nK-i  ^  ~ 


K-1 


} 


ie,  that  the  probability  of  is  independent  of  ipK  prior  to  any  hypothesis  being 

made  about  the  target  and  propagation  path  associated  with  -ipx,  the  following  is  obtained 


'ni-.nif 


ni..nK-i 


K-1 


ie. 


P{A 


mi  ..TUK 

n\..nK 


p{yT 


yT 


y^K-: 


.mi..mK-i  \  orriK  p  /  5,mi..mK--i 
Ani..nK--i  )  Pk  ^  \'^ni..nK-i 


xiffi-l  j 


•®'tl-riK-l  Mk  -  -  - 

E  E 

71^  =  1  mK-  =  l 


TTl\ 


rriK-i 

iJ/K-1 


\  /orriK 

Ani-riK-i  J  Pk 


may  be  possible  to  avoid  this  assumption,  which  has  been  found  to  introduce  an  ordering  property 
on  the  hypothesis  probabilities,  that  is,  the  magnitude  of  the  probabilities  is  afftTcted  by  the  order  in  which 
the  tracks  are  introduced  to  the  algorithm.  A  new  derivation,  which  aims  to  remove  the  assumption  is  in 
progress  and  will  be  presented  in  a  subsequent  publication. 
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If  it  is  assumed  that  no  two  resolved  ground  tracks  which  are  due  to  the  same  target  can 
be  associated  with  the  same  propagation  path  then 


p  j.  ^  1^1  )  Pk  ^  /  ^2) 


Bni..nK_i  Mk  /  - 

E  E  P  ivT 

nK:=l  mK-=l, 


i/1  >")i//C'-i  1 ‘^i..nK-inK 


)0. 


(OlK 

K 


rriK^SK 


Now,  let  us  consider  how  the  likelihood  p  >  -i  Vk-i'^  i  )  is 

calculated.  If  uk  nj,  j  —  —  1  ie,  iI^k  represents  a  new  target,  we  have  no  prior 

information  regarding  what  the  state  of  should  be.  Hence  the  likelihood,  is 

given  by: 


A5« 


PyyK  rK  ’Pif  >Pl  ’"jP/C-I  )''ni..nff_i  ) 


(3) 


V, 

where  Vs  is  the  “volume'^  in  state  space  in  which  the  target  may  be.  The  above  equation 
is  equivalent  to  saying  that  the  target  is  equally  likely  to  be  anywhere  in  the  volume.  Note 
that  this  “volume”  represents  all  the  possible  values  of  the  state  vector  [r  r  a  d]^  where 
r  is  the  range  of  the  target  in  ground  coordinates,  and  a  is  the  azimuth  of  the  target  in 
ground  coordinates. 


Now  consider  the  case  where  tik  =  rij  for  one  or  more  of  j  =  1,  ..,i^  —  1  ie,  ipK 
represents  the  same  target  as  at  least  one  previously  hypothesized  target.  Let  us  assume 
that  the  hypothesis  is  associated  with  T  targets  where  1  <  T  <  K—1  .  Assume 

now  that  the  target  that  tk  corresponds  to  is  tj,  1  <  i  <  T  and  that  prior  to  considering 
there  have  been  h  track  estimates  associated  with  the  current  target  for  this  time 
instant.  Let  the  previous  estimate  of  the  state  of  target  tj  be  and  its  covariance  P-^. 
Since  it  is  assumed  that  (with  covariance  )  is  from  target  ti ,  the  best  available 
estimate  of  y^^  ,  prior  to  actually  determining  y^^  from  tk  and  the  assumed  propagation 
path  mK ,  is  obtained  by  noting  that 


•  all  process  and  measurement  noises  are  assumed  Gaussian, 

•  y^^  is  computed  using  a  Kalman  filter  (which  is  an  unbiased  estimator  if  the  process 
and  measurement  models  are  correct)  and  a  subsequent  ionospheric  transformation 
which  is  assumed  unbiased, 

•  is  computed  using  an  unbiased  estimator, 

•  y^^  and  x^  are  both  estimates  of  the  same  target  at  the  same  time  instant. 


Assuming  all  the  above,  the  best  estimate  (ie,  the  expected  value)  of  y^^  is 

^  xVk  rK  'Pi  f’i/K'-l  I ''ni-riK’-i  /  “ 

The  covariance  of  (y^^  —  x^)  is 


pmK  A  p 


mK 


Vk 


rriK  ~h 


Vk 


]'} 


{Vk^^  -  Xi)  -  {x^  -  Xi)  (y^'^  -  Xi)  -  (x^  -  Xi) 
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Let  =  Xi  —  x^  and  =  Xi  —  ,  where  Xj  is  the  true  state  of  target  ti ,  then 


T^K  =  E 


Xi 


“rriK 

Vk 


i’t-vT 


=  E  |xf  (xj)'  -  xf  ®»)'  -  ij"  (x?)'  +  (sj")'} 

=  £|x!“(x‘)'}-£{x'‘(ir)'}--®{»«''  (*?)'} 

Now.  Pt^E  {x?  (x?)'}  ,  p;?"  4  £  {jij*  ,  pI^'<  4  E  {x?  (yj")'}  .  A 

E  1  hence 


_  ph  _  phtTiK 
K  ^  t  iK 


K 


(4) 


Then  the  likelihood 


where 


Ar  =  v(vT\«T,*l^.yT'<--,vTA'X;:Z‘‘s:) 

=  Af(y“«;xf.7^«) 

=  |2x7^''|-jexp{-i(rr)'(ir'')’‘‘'r} 


(5) 


As  mentioned  earlier,  there  are  some  errors  in  some  of  the  static  equations  presented 
in  [17]  and  [19].  The  key  equations  in  question  are  the  counterparts  of  equations  2,  3,  and 
5  in  this  report.  For  a  description  of  these  errors,  the  reader  is  referred  to  Appendix  A. 

Consider  now  the  estimate  x^'^^  of  Xj  and  its  covariance  obtained  by  the  fusion 
of  x^^,  h  =  1,2, H  —  I  with  y'^^  .  To  derive  this  we  shall  use  the  fundamental  equations 
of  linear  estimation  as  shown  on  pages  44  and  125  of  [4] ,  ie, 

X  =  E  {x\z]  =  X  +  {z  -  z)  (6) 

Exx\z  ~  E  {(x  —  x)  (x  —  x)  —  Pxx  ~  PxzPzz  E zx 

where  x  is  the  random  vector  to  be  estimated,  z  is  the  new  measurement  or  the  observation, 
X  is  the  previous  estimate  of  x,  z  is  the  previous  estimate/prediction  of  z,  x  =  F[x|z]  is 
the  conditional  mean  of  x  given  z,  and 

=  F  [(x  —  x)  (x  —  x)^] 

=  E  [(x  —  x)  (z  —  z)^] 

=  £[(2-x)(x-i)'] 

4  P[(z-f)(z-f)'] 

=  F  [(x  —  x)  (x  -  x)^  |z] 


zx 

Px.- 
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Now,  following  the  approach  shown  in  sections  8.3.3  and  8.4.4  of  [5],  we  replace  the 
terms  in  the  fundamental  equations  of  linear  estimation  with  their  equivalent  terms  in  our 
problem  as  follows 


X 

X 

z 

z 


./l+l 

i 

h 

i 


vT 


[ 


Xi  -  a;/ 


ft.  ^  £  I  [vT  -  *?]  [*i  -  *?]  ’ } 


4]  [«"-*?]'}  “A" 


phniK  pmKk  I  pmji 


K 


Expanding  the  right  hand  side  of  the  replacement  for  Pxz  gives 


{ h  -  4]  [vn 


E 


I  li  -  Ij"]  -  li  -  if  + 1( 


Hence  the  replacement  for  Pxz  is 


ph  _  phmK 


Pxz  -^Pt-  PiK^" 


Expanding  the  right  hand  side  of  the  replacement  for  Pzx  gives 


^i\  \Xi-Xi 


}  =  ®  I  Vk"  -  li  -  I?  +  I( 

=  E I  (if -y^)  (*?)'} 

=  £  {if  (if)' -!/r  (*?)'} 


/'I 

-h 

1 

Xi  - 

J 

and  hence  the  replacement  for  Pzx  is 

Pz: 


_  ph  pmiih 

—  M  ^Ki 


ph  prriKh 

M  ^Ki 
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The  resulting  fusion  equations  are  hence 

xf+i  =  i?  +  [Pt  -  Ptr)  (17")-'  {vT  - ^?) 

pM  ^pk_  (^pk  _  pk^K^  (^pk  _  p<nKk-\^ 

Note  that  the  cross-covariance  matrices  and  are  extremely  difficult,  if 

not  impossible,  to  determine.  It  seems  that  it  may  be  more  appropriate  to  use  the  actual 
measiurements  from  which  the  (radar  coordinate)  tracks  are  derived,  than  attempt  to 
rigorously  compute  the  cross-covariances  of  the  track  estimates.  Using  measurements 
would  remove  some  of  the  contributions  to  the  cross-covariance  terms  (this  will  be  discussed 
later  in  the  report).  As  a  consequence,  the  approach  that  is  taken  here  is  to  approximate 
the  values  of  the  cross-covariance  matrices.  Later,  an  algorithm  for  the  creation  of  the 
hypothesis  tree  using  the  measurements  will  be  derived,  the  aim  being  to  compare  the 
comparative  performance  of  the  two  approaches. 

As  described  above,  the  static  MPTF  algorithm  computes  the  association  probabilities 
and  performs  fusion  for  every  possible  combination  of  target  and  propagation  path  associ¬ 
ation.  Unfortunately,  for  the  numbers  of  tracks  and  propagation  paths  found  in  practice, 
the  algorithm  cannot  be  implemented  without  some  form  of  hypothesis  pruning.  The  same 
applies  to  the  DMPTF  algorithm  as  it  uses  the  “static”  MPTF  algorithm  for  its  initial 
creation  of  its  hypothesis  tree.  An  efficient  algorithm  for  performing  recursive  pruning  of 
the  MPTF  hypothesis  tree  has  been  developed  and  is  described  in  section  4,  following  a 
discussion  regarding  the  computational  complexity  of  MPTF  in  section  3. 


(7) 

(8) 
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3  Computational  Complexity  of  Multipath 

Track  Fusion 


In  the  absence  of  a  hypothesis  management  strategy,  the  number  of  association  hy¬ 
potheses  produced  by  the  MPTF  algorithm  can  quickly  become  too  large  for  real  time 
execution.  Some  preliminary  work  on  overcoming  this  difficulty  was  described  in  Refs. 
[18]  and  [19],  but  the  approach  presented  was  found  to  be  of  very  limited  usefulness.  This 
section  presents  equations  for  the  number  of  hypotheses  that  the  MPTF  algorithm  gener¬ 
ates  and  then  gives  a  summary  of  the  pruning  work  described  in  [18]  and  [19],  followed  by 
reasons  for  the  inadequacy  of  the  approach. 

Consider  the  path-independent  hypothesis  tree  for  a  cluster  of  J  tracks.  The  number 
of  track-to-target  association  hypotheses  at  a  given  depth  J  of  the  hypothesis  tree  for 
the  path  independent  case  was  evaluated  in  [17]  to  be  the  Bell  number  (or  exponential 
number)  [6] 

J 

T=1 

where  S  {J,  T)  is  the  Stirling  number  of  the  second  kind  and  is  given  by 

=  i<T<j 

■  t=l 


where 


T\ 

(T-t)W.' 


Now  consider  the  propagation  paths  that  are  to  be  associated  with  radar  tracks.  If  all 
combinations  of  ionospheric  paths  are  allowed,  it  can  be  easily  deduced  that  the  number 
of  association  hypotheses  at  a  given  depth  J  of  the  path  dependent  hypothesis  tree,  with 
M  available  ionospheric  paths,  is 

J 

r=i 


As  an  example  of  what  can  be  expected  in  real  surveillance  scenarios,  for  a  cluster 
of  10  tracks  and  9  possible  ionospheric  paths  the  number  of  path-dependent  hypotheses 
in  the  final  (i.e.,  10th)  level  of  the  hypothesis  tree  is  (10, 9)  «  4.04  x  10^^.  The  total 
number  of  hypotheses  that  would  need  to  be  evaluated  to  build  the  tree  is 

J 


For  10  tracks  and  9  possible  ionospheric  paths  H{J,M)  w  4.  13  X  10^^.  Clearly  the 
hypotheses  cannot  be  exhaustively  evaluated  in  real  time  nor  stored  in  computer  memory. 
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An  approach  to  reduce  the  number  of  association  hypotheses  a  priori  was  described  in 
Refs.  [18]  and  [19].  In  that  scheme,  the  formulation  of  hypotheses  relied  on  the  application 
of  physical  (deterministic)  constraints. 

The  assumption  that  the  physical  constraints  are  based  on  is  that  two  radar  tracks 
having  different  range  estimates  may  arise  from  the  same  target  only  if  the  track  at  great¬ 
est  range  is  associated  with  the  propagation  path  of  greatest  path  length.  This  prior 
information  can  be  used  for  the  a  priori  pruning  of  hypotheses. 

The  physical  constraint  is  apphed  as  follows.  Tracks  tj  ,  j  =  1, ..,  J  are  ordered  with  re¬ 
spect  to  one  of  the  components  of  the  track  state  vector  xpj ,  typically  the  range  component 
Rj  which  gives 

Ri  <  R2  <  ...  <  Rj.  (9) 


Likewise,  the  propagation  paths  ruj  ,  for  the  jth  track  j  =  1, ..,  J  are  ordered  with  respect 
to  resulting  ground  range  rj'^  so  that 


m] 

r/ 


M,- 


<  r,-  <  ...  <r. 


(10) 


For  path-dependent  hypotheses  of  the  form  where  rii  =  nj  ,  i  <  j  ,  the  ordering 

of  the  tracks  and  propagation  paths  in  conjunction  with  the  above  assumption  leads  to 
the  constraint 


rrii  <  rUj.  (11) 

This  approach  promises  the  advantage  of  decreeising  the  number  of  hypotheses  that  must  be 
evaluated,  hence  minimizing  computational  expenditure.  Unfortunately,  as  demonstrated 
below,  at  best,  this  is  only  a  partial  solution,  as  it  does  not  prune  enough  hypotheses  to 
make  the  MPTF  algorithm  practical  for  OTHR  applications. 

When  the  deterministic  constraint  is  applied  with  J  tracks  and  with  M  available 
ionospheric  propagation  paths,  the  number  of  path-dependent  association  hypotheses  is 
given  by 

J-l 

H{J,M)=  C'/-iCf_j^(j,M), 

j=max(0,  J— A/) 

where  the  recursive  expression  is  initialised  with  H{0,  M)  =  1.  A  table  giving  values 
of  H{J,M)  evaluated  for  1  <  J  <  6  and  1  <  M  <  8  has  been  presented  in  Ref.  [18]. 
Returning  to  the  earlier  example,  with  the  number  of  tracks  J  —  10,  the  number  of 
ionospheric  propagation  paths  M  =  9,  and  the  constraint  applied,  the  final  level  of  the 
hypothesis  tree  has  H{10, 9)  =  4.89  x  10^^  hypotheses.  The  total  number  of  hypotheses  is 
H  =  5.08  X  10^^  ,  which  is  approximately  100  times  less  than  the  number  of  hypotheses 
in  the  absence  of  any  a  priori  pruning.  However,  it  is  clear  that  the  application  of  the 
deterministic  constraint  alone  does  not  reduce  the  number  of  hypotheses  anywhere  near 
enough  to  make  the  MPTF  algorithm  practical.  Additional  hypothesis  pruning  methods 
are  required  regardless  of  whether  deterministic  constraints  are  applied.  There  are  also 
other  difficulties  with  using  physical  constraints  which  make  the  approach  impractical  even 
for  preliminary  pruning.  These  are  described  below. 
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Because  the  ordering  of  tracks  and  propagation  paths  is  performed  using  estimates  of 
random  variables,  applying  the  constraint  (11)  can  give  rise  to  errors.  Firstly,  for  two 
associated  radar  tracks  arising  from  propagation  via  the  “mixed  paths”  such  as  EF  and 
FE,  the  track  range  is  nominally  identical.  Thus  the  range  orderings  represented  by  (9) 
and  (10)  may  result  in  a  correct  path-dependent  hypothesis  not  being  formulated  when 
(11)  is  applied.  A  simple  method  to  overcome  this  difficulty  may  be  formulated  whereby  a 
parameter  6R  which  is  a  function  of  the  track  range  variance  is  used  to  apply  the  constraint 
(11)  when 


Ri  +  SRi  <  Rj  —  5Rj  i<j.  (12) 

However,  this  reduces  the  already  limited  amount  of  pruning  that  is  achieved. 

Secondly,  the  software  implementation  of  the  physical  constraint  is  found  to  be  com¬ 
plicated  when  the  coordinate  registration  advice  is  considered,  eg., 

•  The  CR  system  may  advise  a  different  set  of  propagation  paths  for  different  tracks  in 
which  case  the  ordering  of  paths  according  to  a  given  track  may  not  be  appropriate 
to  the  paths  advised  for  another  track.  Note  that  in  (9)  the  propagation  paths  are 
ordered  according  to  the  transformations  for  the  track  only. 

•  When  the  number  of  tracks  and  propagation  paths  and  their  order  changes  with  time, 
the  changing  order  must  be  somehow  accommodated  in  formulating  the  updated 
hypotheses. 

This  introduces  unwarranted  complexities  in  computation  and  management  of  data  as  well 
as  coding  complexities.  It  is  important  to  note  that  the  original  aim  of  a  priori  hypothesis 
pruning  is  to  reduce  computation  and  data  management  when  formulating  and  evaluating 
hypotheses. 

Finally,  from  a  mathematical  standpoint  the  physical  constraint  approach  is  essentially 
a  deterministic  tack-on  to  a  stochastic  algorithm  (MPTF).  This  may  have  been  justifiable 
had  it  achieved  substantial  computational  efficiencies;  however,  this  was  found  to  not  be 
the  case. 

For  these  reasons  it  was  decided  to  abandon  the  physical  constraint  approach  sum¬ 
marized  above  in  favour  of  the  recursive  a  posteriori  pruning  algorithm  described  in  the 
following  sections.  Note  that  a  simpler  constraint,  rrij  ^  mj  if  Ui  =  Uj,  i  ^  j  has  been 
retained  where  appropriate.  In  words,  the  constraint  states  that  two  tracks  due  to  the 
same  target  cannot  be  propagating  via  the  same  propagation  path.  This  constraint  does 
not  suffer  from  the  disadvantages  of  (11). 
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4  Recursive  Hypothesis  Pruning 

In  order  to  maJce  MPTF  computationally  feasible,  a  new  pruning  algorithm  has  been 
developed  which  exploits  the  recursive  nature  of  the  MPTF  algorithm  by  pruning  unlikely 
hypotheses  at  each  level  of  the  hypothesis  tree  as  it  is  built,  prior  to  creation  of  the  next 
level.  The  decision  as  to  which  hypotheses  to  prune  is  made  on  the  basis  of  hypothesis 
probabilities  which  are  calculated  using  all  components  of  the  track  estimate  state  vectors 
and  path  transformations.  If  H  is  the  number  of  hypotheses  formulated  by  the  MPTF 
algorithm  at  a  particular  level  of  the  hypothesis  tree,  then  the  pruning  scheme  retains  the 
K  best  (highest  probability)  hypotheses.  These  hypotheses  are  subsequently  used  in  the 
MPTF  algorithm  to  form  the  next  level  of  the  hypothesis  tree.  The  remaining  H  —  K 
hypotheses  are  pruned  and  therefore  are  not  used  in  forming  hypotheses  in  higher  levels 
of  the  tree.  The  integration  of  the  pruning  scheme  with  the  MPTF  algorithm  is  outlined 
in  the  following  pseudo-code,  with  the  last  two  steps  being  the  additional  steps  required 
to  perform  hypothesis  pruning. 

-  Using  the  first  track  estimate,  formulate  all  path- dependent  hypotheses  and  calculate 
their  probabilities; 

-  While  unprocessed  track  estimates  remain: 

-  Select  the  next  track  estimate; 

-  Formulate  new  path- dependent  hypotheses  using  the  current  track  estimate  and  the 

path- dependent  hypotheses  remaining  in  the  previous  level  of  the  hypothesis  tree; 

-  Calculate  probabilities  of  all  new  path- dependent  hypotheses; 

-  Determine  and  mark  the  K  most  probable  hypotheses; 

-  Remove  all  unmarked  hypotheses  from  the  current  level  of  the  hypothesis  tree; 

The  major  advantages  of  the  new  approach  are: 

1.  It  enables  a  huge  reduction  in  the  number  of  hypotheses  that  are  evaluated,  and  easy 
tailoring  to  the  computational  capabilities  of  the  computer  on  which  the  algorithm 
is  to  be  implemented  by  adjusting  the  value  of  K. 

2.  The  approach  is  consistent  with  the  probabilistic  model  being  used  to  model  the 
system,  using  all  the  components  of  the  state  vector  and  path  transformations,  and 
their  covariances. 

The  approach  has  a  potential  disadvantage,  i.e.,  finding  the  K  best  out  of  an  initial 
H  hypotheses  is  computationally  expensive  when  both  K  and  H  are  large  if  a  simple 
sequential  search  algorithm  is  used.  In  order  to  overcome  this  potential  limitation,  several 
algorithms  for  finding  the  K  best  hypotheses  were  investigated  with  the  aim  of  producing 
a  computationally  efficient  pruning  algorithm.  The  following  two  subsections  give  an 
outline  of  the  algorithms  that  were  investigated  and  the  results  of  comparative  testing, 
respectively. 


4.1  K-best  Search  Algorithms 
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Prior  to  describing  the  algorithms  that  were  investigated,  it  is  appropriate  to  make  a 
distinction  between  hypothesis  management  in  multiple-hypothesis  tracking  (MHT)  [23] 
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and  that  required  in  MPTF.  Techniques  such  as  Lagrangian  Relaxation  [20],  which  have 
been  applied  for  finding  the  best  set  of  association  hypotheses  in  MHT,  are  not  applicable 
for  MPTF.  This  is  because  MPTF  is  not  cast  in  the  form  of  an  assignment  problem,  as 
is  the  case  for  MHT.  In  the  assignment  problem  there  is  competition  for  “resources”;  for 
example,  in  MHT,  tracks  compete  for  measurements.  In  the  case  of  MPTF,  there  is  no 
such  competition.  As  a  result,  an  alternative  algorithm,  to  essentially  solve  a  different 
problem,  needed  to  be  found.  The  development  of  this  algorithm  is  outlined  below. 

Data  structures  and  algorithms  from  the  computer  science  literatiure  were  investigated 
to  yield  schemes  that  would  efficiently  find  the  K  hypotheses  with  highest  probability 
from  a  set  of  if  initial  hypotheses  (ie,  perform  a  K-best  Search).  Several  of  the  most 
promising  techniques  were  examined  in  detail  and  compared,  with  the  aim  of  producing  a 
computationally  efficient  pruning  scheme  to  be  implemented  within  the  MPTF  algorithm. 
The  algorithms  and  data  structures  that  were  investigated  and  tested  in  detail  were  FIND 
[9],  AVL  Tfrees  [1],  [2],  Skip  Lists  [21]  and  Binary  Heaps  [26],  [10].  FIND  is  generally  used 
in  sorting,  whereas  AVL  Ttees  and  Skip  Lists  are  generally  used  for  searching  algorithms, 
and  the  Binary  Heap  for  constructing  Priority  Queues  [26],  [10].  In  addition  to  the  above 
techniques,  the  performance  of  a  repeated  Sequential  Search  was  measured.  This  was  used 
as  a  reference  against  which  the  more  sophisticated  approaches  could  be  compared. 

The  algorithms  that  were  examined  in  detail  are  outlined  in  the  following  paragraphs. 
In  the  algorithm  descriptions,  conventional  computer  science  terminology  is  employed 
when  appropriate.  Specifically,  a  record  is  defined  as  an  object  which  stores  information. 
A  key  is  a  component  of  a  record  which  can  be  used  to  reference  the  record.  In  a  structure 
that  has  many  records,  their  keys  can  be  used  to  find  individual  records  (ie  search)  or  to 
order  the  records  (ie  sort).  In  the  MPTF  pruning  algorithm  the  records  correspond  to  the 
hypotheses  and  the  keys  are  their  probabilities. 

Sequential  Search  The  simplest  algorithm  for  finding  the  best  K  of  H  records  is 
to  perform  a  Sequential  Search  K  times.  At  each  pass,  the  best  record  is  found  and  then 
tagged  so  that  it  is  not  considered  in  subsequent  searches.  The  algorithm  is  simple  but 
inefficient  when  H  and  K  are  large,  having  a  computational  complexity  of  0{KH).  A 
summary  of  the  algorithm  follows: 

Starting  with  a  list  of  H  records,  repeat  the  following  K  times: 

-  Scan  through  the  list  of  records  to  find  the  record  with  greatest  key; 

-  Record  and  tag  that  record  so  it  is  not  scanned  in  subsequent  passes; 

FIND  Sorting  techniques  based  upon  Hoare’s  FIND  algorithm  [9]  are  widely  used, 
for  example  in  finding  the  median  of  a  set  of  numbers.  Consider  an  array  A(1),A(2), 
...  ,A{K),  ...  ,A{H)  containing  the  keys  of  H  records.  FIND  partitions  the  data  such 
that  A(l) . . .  A{K  —  1)  >  A{K)  and  A{K  +  1) . . .  A{H)  <  A{K),  but  does  not  require 
either  partition  to  be  fully  ordered. 

AVL  Tree  There  are  many  binary  tree  algorithms  which  seek  to  ensure  that  an 
efficient  search  tree  is  always  maintained  by  balancing  the  tree  [26].  An  example  of  a 
balanced-tree  algorithm  is  the  AVL  tree  developed  by  Adelson-Velskii  and  Landis  [1],  [2]. 
The  AVL  tree  is  a  binary  tree  which  has  the  property  that  the  sub-trees  of  every  node 
differ  in  height  by  at  most  one.  A  balance  factor  is  added  to  each  node  in  the  tree  and 
rotations  of  nodes  within  the  tree,  based  on  the  balancing  factor,  are  performed  to  retain  a 
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balanced  tree.  While  balanced-tree  algorithms  ensure  good  performance,  they  are  complex 
to  code  and  the  overheads  in  maintaining  a  balanced  tree  are  significant. 

Skip  List  A  Skip  List  is  an  ordered  list  data  structure  with  additional  forward  pointers 
to  skip  multiple  records  in  the  list,  thereby  reducing  the  time  taken  to  traverse  the  list  [21]. 
Each  record  is  represented  by  a  node  in  the  skip  list.  A  node  has  level  i  if  it  has  i  forward 
pointers.  A  random  number  generator  is  used  to  specify  the  level  of  each  node  according 
to  the  skip  list  parameter  p,  which  is  the  probability  that  a  node  has  a  level  i  pointer 
given  that  it  has  a  level  i  —  1  pointer.  The  parameter  p  is  chosen  to  maximise  algorithm 
performance,  often  set  to  1/4  or  1/2.  The  average  performance  of  the  Skip  List  is  in  the 
order  of  KlnH,  but  its  performance  reduces  to  that  of  the  repeated  Sequential  Search 
in  the  worst  case.  The  poor  worst  case  performance  is  considered  a  serious  limitation  for 
real-time  applications. 

Binary  Heap  A  Binary  Heap  [26]  is  a  special  type  of  binary  tree,  represented  as  an 
array,  which  maintains  the  first  element  in  the  array  as  the  record  with  the  minimum  (or 
maximum)  key.  The  structure  is  arranged  such  that  the  children  of  node  j  in  the  array 
are  in  position  2j  and  2j  -f  1,  and  the  parent  of  the  jth  node  is  in  position  j  div  2.  The 
elements  in  the  structure  satisfy  the  Heap  property  ie,  the  children  of  node  j  have  larger 
(or  smaller)  keys  than  the  record  at  node  j.  The  Binary  Heap  is  often  used  as  a  Priority 
Queue  algorithm  [10],  [26]  because  it  allows  very  simple  and  efficient  implementation  of 
the  Priority  Queue  insertion,  deletion  and  replacement  operations. 

In  this  work,  three  Priority  Queue  operations,  insert,  delete  and  replace,  are  supported 
by  the  Binary  Heap  in  order  to  implement  the  AT-best  Search  algorithm.  The  insert 
operation  places  a  new  record  into  the  Heap  after  the  previous  last  record,  and  then 
reorders  the  records  so  that  they  again  satisfy  the  Heap  property.  The  delete  operation 
removes  the  top-most  (i.e.  the  minimum  key  and  first  index  in  the  array)  record  from  the 
heap,  and  then  reorders  the  remaining  records  as  necessary.  The  replace  operation  replaces 
the  top-most  record  in  the  heap  with  the  new  record  and  then  reorders  the  remaining 
records.  A  detailed  description  of  the  three  operations  is  given  in  Ref.  [26]. 

The  repeated  Sequential  Search  and  FIND  were  implemented  by  entering  all  H  records 
into  a  one-dimensional  array  and  then  determining  the  highest  K  records  (ie,  the  K 
records  with  the  highest  keys).  With  AVL  Trees,  Skip  Lists  and  Binary  Heaps  there  were 
two  options  available,  owing  to  the  data  structures  being  used.  Because  the  three  data 
structures  incorporate  an  ordering  mechanism  by  virtue  of  their  construction,  the  act  of 
inserting  a  record  into  each  structure  ensures  that  records  within  the  structure  are  ordered 
by  their  keys.  AVL  Trees  and  Skip  Lists  are  fully  ordered  structures  whereas  Binary  Heaps 
are  only  partially  ordered.  Most  importantly,  the  minimum  record  in  the  structure  can 
easily  be  accessed  and  replaced.  Hence,  for  these  three  cases,  two  types  of  A-best  Search 
algorithm  were  implemented.  The  first  type  sequentially  inserts  all  H  records  into  the 
structure  (as  was  done  for  Sequential  Search  and  FIND),  then  extracts  the  highest  K 
records  from  the  structure  one-by-one.  The  second  approach  stores  only  K  records  in 
the  structure  at  any  one  time,  by  monitoring  the  minimum  record  in  the  structure  and 
only  inserting  those  records  that  have  a  key  which  exceeds  that  minimum.  The  second 
algorithm  is  summarised  by  the  following  pseudo-code: 
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For  the  first  K  records  in  the  list  of  H  records: 

-  Add  the  records  into  the  data  structure. 

For  all  subsequent  records: 

If  the  key  of  the  selected  record  exceeds  that  of  the  minimum  record  already  in  the 

structure,  then 

-  Remove  the  minimum  key  record  from  the  structure. 

-  Add  the  new  record  into  the  structure. 

A  key  point  to  note  about  the  second  approach  is  that  it  reduces  the  number  of 
insertions  into  the  structure  as  well  as  minimizing  memory  requirements. 

4.2  Performance  Comparison 

The  average  performance  of  Sequential  Search,  FIND,  AVL  IVee,  Skip  List  and  Binary 
Heap  schemes,  when  used  as  AT-best  Search  algorithms,  was  compared  by  implementation 
(in  C++)  and  testing  on  a  233  MHz  Pentium  PC,  running  Windows  NT  4.0.  Input  data 
was  simulated  by  generating  1000  sets  of  H  random  numbers,  with  the  aim  of  determining 
the  K  highest  numbers  in  each  set.  The  measure  of  performance  was  the  average  run  time 
for  executing  the  algorithms,  which  was  determined  as  a  function  of  K  and  H. 

Figure  5  presents  the  average  times  taken  to  find  the  K  highest  values  from  a  set  of  Af 
numbers.  Each  set  of  numbers  was  randomly  generated  and  uniformly  distributed  in  the 
interval  (0  —  1).  The  three  algorithms  based  on  ordered  data  structures,  (i)  Binary  Heap, 
(a)  AVL  tree  and  (Hi)  Skip  List,  were  each  implemented  with  structure  sizes  of  both  H 
and  K  records.  In  the  figure,  as  well  as  the  following  text,  the  structure  size  is  written  in 
parentheses  following  the  structure  name. 

Figure  5(a)  and  (c)  show  that  the  performances  of  the  AVL  Tree  (H)  and  Skip  List  (H) 
algorithms  are  similar  to  that  of  Sequential  search  when  K  is  small  but  are  substantially 
better  for  K  =  500.  The  fastest  algorithms  for  most  values  of  K  and  H  are  FIND,  Binary 
Heap  (K),  AVL  Tree  (K),  and  Skip  List  (K).  Consistently  high  performance  is  attained 
by  the  FIND  and  Binary  Heap  (K)  algorithms,  especially  for  K  —  500,  making  them  the 
best  candidates  for  use  in  the  hypothesis  pruning  algorithm.  Of  the  two.  Binary  Heap 
(K)  is  considered  a  better  choice  because  it  requires  an  array  of  size  K,  in  contrast  to 
the  FIND  algorithm  which  requires  an  array  of  size  H.  It  should  also  be  noted  that  the 
Binary  Heap  (K)  algorithm  is  simple  to  code,  which  is  an  important  consideration  in  the 
implementation  of  the  multipath  track  fusion  algorithm. 

For  the  reasons  given  above,  the  Binary  Heap  (K)  algorithm  was  selected  for  the 
implementation  of  hypothesis  pruning  in  MPTF. 
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*  Sequential  Search  •  Find 

♦  Binary  Heap  (H)  ■  Binary  Heap  (K) 

+  AVL(H)  *AVL(K) 

▼  Skip  List  (H)  ^  Skip  List  (K) 

Figure  5:  Comparison  of  average  run-time  of  candidate  K-best  Search  schemes  as  a  func¬ 
tion  of  the  initial  number  of  records,  H,  for  the  cases  when  K  =  50,  K  —  100  and 
K  —  500.  Plots  (a),  (c),  and  (e)  show  the  same  data  as  (b),  (d)  and  (f),  respectively,  but 
with  different  scales  on  the  vertical  axes. 
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5  Hypothesis  Pruning  Algorithm 
Implementation 

As  stated  in  the  previous  section  the  pruning  algorithm  that  was  found  to  be  most 
efficient  was  the  one  based  on  a  Binary  Heap  structure.  This  algorithm  was  coded  on  the 
MPTF  testbed  and  its  implementation  was  found  to  extend  the  capability  of  MPTF  to 
much  larger  numbers  of  tracks  and  propagation  paths  than  previously  possible,  making 
it  feasible  to  use  MPTF  without  further  extensions  when  there  is  a  moderate  number 
of  tracks  and  propagation  paths,  ie,  more  than  20  tracks  with  10  propagation  paths  with 
computational  power  presently  available  to  the  MPTF  algorithm  in  the  JFAS  radar.  While 
the  OTHR  can  be  expected  to  have  considerably  more  than  20  tracks  over  its  entire 
coverage  area,  only  tracks  which  are  relatively  close  to  one  another  can  conceivably  be  due 
to  the  same  target,  the  multiple  tracks  (from  a  single  target)  forming  clusters  of  almost 
invariably  less  than  9  tracks  and  usually  4  tracks  or  less.  This  characteristic  was  used  to 
enable  MPTF  to  deal  with  very  large  numbers  of  tracks  (eg.,  hundreds  of  tracks  with  ten 
propagation  paths).  This  was  achieved  by  the  development  of  a  clustering  algorithm  which 
groups  tracks  which  can  conceivably  be  due  to  the  same  target  and  then  passes  on  these 
groups  of  tracks  to  the  MPTF  algorithm  for  processing.  The  purpose  of  the  clustering 
algorithm  is  not  to  decide  which  tracks  are  due  to  the  same  target,  but  instead,  to  form 
groups  of  tracks  which  may  be  due  to  the  same  target,  while  maintaining  group  size  to  less 
than  about  say  20  tracks,  and  at  the  same  time  not  placing  tracks  from  the  same  target 
into  separate  clusters.  The  combination  of  pruning  and  clustering  was  implemented  on 
the  testbed  and  also  a  prototype  version  was  incorporated  in  the  actual  Jindalee  OTHR. 
The  combination  worked  very  well,  handling  tracks  in  the  entire  coverage  area  of  the  radar 
effectively  and  without  creating  unacceptable  delays  in  computation.  No  further  details 
of  the  clustering  algorithm  will  be  given  here  as  it  will  be  the  subject  of  a  future  report. 

As  mentioned  in  an  earlier  section,  the  pruning  technique  that  has  been  developed 
exploits  the  recursive  nature  of  the  MPTF  algorithm  by  pruning  low  probability  hypotheses 
at  each  recursion,  thus  avoiding  the  exponential  growth  in  the  number  of  hypotheses. 
A  possible  disadvantage  of  this  approach  is  that  hypotheses  that  are  unlikely  in  early 
recursions  may,  under  some  circumstances,  subsequently  lead  to  the  creation  of  likely 
hypotheses,  had  pruning  not  been  performed.  Preliminary  investigations  have  indicated 
that  this  is  not  likely  to  be  a  problem  and  that  the  MPTF  algorithm  will  perform  very 
effectively  with  quite  heavy  pruning  implemented  in  the  manner  described.  While  no 
detailed  studies  have  yet  been  performed  to  determine  the  best  choice  for  the  number  of 
hypotheses  that  should  be  retained,  in  the  investigations  to  date,  retention  of  the  best  100 
hypotheses  in  each  algorithm  recursion  has  proven  quite  successful. 

Finally,  it  should  be  noted  that  the  hypothesis  pruning  algorithm  that  has  been  de¬ 
veloped  here  for  OTHR  multipath  track  fusion  may  also  be  applied  to  multisensor  track 
fusion  algorithms  that  have  been  cast  in  the  multihypothesis  framework. 


23 


DSTO-RR-0223 


6  Dynamic  Multipath  Track  Fusion  Algorithm 

In  earlier  work,  as  described  in  [17]  and  [19],  the  procedure  for  update  A:  =  0  is 
simply  repeated  for  all  subsequent  updates  {k  =  1, 2, 3, ...),  with  each  update  being  treated 
independently.  The  major  disadvantages  of  this  approach  are 

•  only  synchronous  estimates  can  be  included  in  the  hypothesis  tree,  hence  multisensor 
tracks  cannot  be  accommodated, 

•  of  more  fundamental  theoretical  importance,  temporal  behaviour  of  tracks  is  com¬ 
pletely  ignored, 

•  when  multipath  tracks  drop  out  or  come  into  existence  there  is  no  theoretically 
justifiable  way  to  “link”  hypotheses  from  consecutive  updates. 

In  [19],  possible  extensions  to  include  target  dynamics  are  outlined.  It  is  suggested  that 
the  tracker’s  target  dynamic  model  be  used  to  perform  predictions  (in  radar  coordinates), 
which  would  then  be  transformed  to  ground  coordinates.  The  extensions  only  involve 
synchronous  updating  of  the  tracks  that  are  used  as  inputs  to  the  MPTF  algorithm. 
Of  considerably  more  importance,  there  is  no  description  of  how  hypothesis  probability 
updating  can  be  performed,  other  than  the  presentation  of  two  very  general  equations  that 
have  been  transcribed  from  [14].  These  equations  offer  little  insight  into  how  hypothesis 
probability  updating  would  actually  be  performed.  Initiation  and  termination  of  tracks  is 
also  discussed  but  no  solution  presented. 

An  algorithm  which  overcomes  all  of  the  above-mentioned  limitations  is  derived  and 
presented  in  detail  in  the  following  subsections. 


6.1  Incorporation  of  Target  Dynamics 

The  dynamic  multipath  track  fusion  (DMPTF)  algorithm  which  will  now  be  derived 
overcomes  the  major  limitations  of  the  static  MPTF  algorithm  by  using  an  approach  in 
which  target  behaviour  is  modelled.  This  enables  hypothesis  probability  calculation  and 
fusion  of  estimates  that  are  temporally  separated,  enabling  asynchronicity  to  be  accom¬ 
modated  and  temporal  behaviour  to  be  taken  into  account.  In  the  DMPTF  algorithm, 
the  hypothesis  tree  constructed  at  fc  =  0  is  as  described  earlier  in  section  2.  However, 
instead  of  creating  a  new  tree  at  each  update,  as  is  done  in  the  static  MPTF  algorithm, 
hypothesis  probabilities  and  fused  estimates  of  the  initial  tree  are  updated  as  new  data 
become  available  at  updates  k  >  0. 

The  following  assumptions  are  made  when  developing  the  algorithm: 

1.  If  track  Tj  is  due  to  target  with  propagation  via  path  mj  at  update  k  =  0,  then  it 
remains  so  for  all  future  updates.  This  is  clearly  true  if  the  tracking  algorithm  has 
performed  its  function  correctly. 

2.  The  number  of  tracks  is  constant  for  all  updates.  This  is  only  a  temporary  assump¬ 
tion  that  is  utilized  in  the  initial  development  of  the  DMPTF  algorithm  and  is  later 
removed  by  the  development  of  an  adjunct  algorithm  which  will  be  described  later. 
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As  a  result,  no  new  hypotheses  are  generated  after  the  initial  creation  of  the  tree  at  fc  =  0. 
For  updates  A;  >  0,  only  the  hypothesis  probabilities  and  their  associated  fused  estimates 
change. 

Now  let  us  consider  updates  k  =  1,2, 3, ...  .  The  hypothesis  tree  built  up  at  fc  =  0  can 
be  used  as  the  basis  for  the  tree  at  subsequent  updates  as  described  below. 

In  order  to  shorten  notation,  let  us  define  the  symbol  (k)  such  that  (k)  = 
ipi  (k)  (k).  In  the  creation  of  the  hypothesis  tree  (at  update  A;  =  0)  we  derived 

the  probability  F  (0)  }  of  each  path  dependent  hypothesis  ,  mj  € 

{1, 2, ..,  iWj}  ,  Tij  E  ^1,2, given  all  the  data  available  at  update  ^  =  0 ,  ie 
the  estimates  tpi  (0) , V’J  (0),  as  described  in  section  2.  Now  let  us  define  the  following 
shorthand  notation: 


Dk 


-  (0)  ,..,i)j{0)  ,..,'i/jj{l)  ,....,i{}i{k-l)  ,..,ilij{k-l),tpi  (k)  ,..,il)j  (k) 


ie.,  represents  all  the  (radar  coordinate)  track  data  available  up  to  and  including 
update  k  .  Hence 


p  f  \mi..mj 


D''  }  4  p  |vlr^  (0)  ,  (1)  , ..,  {k)  ] 


Consider  the  updating  of  the  probability  of  hypothesis  when  the  first  track 

estimate  for  update  k  ,  (ie,  (A;))  is  considered.  Using  Bayes  rule  we  have 

*-i  1  _  P  {’Pi  (k)  P*-‘)  P  (Ag-'y  |r>*-‘ } 


p(V'i  {k) 

p  (^1  (fc)  P  } 


Mj 


Ml 


nj=l  hi=lfh.j=l  mi=l 


Now,  as  noted  earlier,  in  the  implementation  of  MPTF  the  ground  co-ordinate  esti¬ 
mates  (k)  are  used  to  perform  the  probabihty  calculations,  hence  we  need  an  equation 
in  terms  of  them.  Hence  again,  wherever  we  encounter  the  conditioning  event  we 
replace  ipj  (k)  with  (k).  Now,  when  we  have  the  event  ^Tii.'.^ir^'^  we  note  that  this 
includes  the  event  (as  well  as  a  number  of  other  conditions).  Hence  the  above  equation 
can  be  replaced  by 


P 


_ p  (vT'  w  p  1^*“* } _ 

Bni..nj_i  1  Mj  Ml 

E  ..  E  E  ■■  E  p{!/f‘WIAf,‘.:,7.B'-‘)p{A?.‘;7iD‘-i} 

nj=l  ni=l?fij=l  ?ni=l 


To  calculate  p  (y^^^  {k)  \XZ\:Z^ ,  we  need  information  regarding  the  dynamics  of 

the  target  that  yj"^  {k)  is  associated  with,  or  predictions  to  be  provided  by  the  individual 
trackers  producing  the  tracks  Tj  .  The  approach  that  will  be  taken  will  be  the  former,  ie, 
using  a  model  of  target  dynamics.  The  reasons  for  choosing  this  approach  are 
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•  Less  information  needs  to  be  provided  to  the  fusion  algorithm  by  the  trackers,  re¬ 
ducing  communication  bandwidth. 

•  Using  target  dynamics  gives  a  more  general  solution,  for  example  if  measurements 
are  fused  instead  of  track  estimates  as  will  be  described  later  in  the  document,  the 
approach  can  still  be  applied. 

Now  let  us  assume  that  target  satisfies  the  following  state  equation: 

Xi{k)  =  F{k-l)xi{k-l)  +  v{k-l)  k  =  1,2,3,....  (13) 

where  F  (k  —  1)  is  the  state  transition  matrix  of  each  of  the  targets  for  the  time  interval 
between  update  A:— 1  and  update  k ,  and  v{k  —  1) ,  A:  =  1, 2, ...  is  the  sequence  of  zero-mean 
white  Gaussian  process  noise  with  covariance  [u  (A:  —  1)  u  (A:  —  1)  ]  =  Qk-i  ■  Hence, 
given  the  fused  estimate  xf  [k  —  1)  for  the  hypothesis  being  considered,  where  H  is  the 
number  of  track  estimates  that  were  combined  at  time  A:  —  1  to  obtain  xf^  (A:  —  1) ,  we  have 
the  following  prediction  at  time  k  based  on  target  dynamics 

x°(A:)  =  F{k-l)xf{k-l)  (14) 

P^{k)  -  F(A:-l)^"(A:-l)F(A:-l)'-t-Qit-i 

Using  the  above  we  can  easily  obtain  the  likelihood 

ArW  = 

=  M{yT‘(k);x?{k),7T'{h)) 

=  \2nTr'  (fc)riexp|-l  {yT‘  (k)  -  xf  (k))’ 7^'  {*)-'  for  (*)-!?(*))} 

where 

7^^  (k)  =  (k)  -  (A:)  -  {k)  +  {k) 

and 

P<‘{k)  4  £{i»  (*)*?«'} 

prik)  ^  Bforwirw'} 

A  £{x»{*)iirw'} 
pr“w  =  E{ff'{k)xf{kY} 

In  the  above 

{k)  =  Xi  {k)  -  x°  {k) 

y-‘(fc)  =  x,{k)-yr{k) 

Consider  now,  the  updating  of  the  probability  of  the  hypothesis  when  the 

track  estimate  for  update  k ,  (ie,  ipj  [k) ,  1  <  j  <  J  )  is  considered.  Using  Bayes  rule 
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we  have 


P  (fe)  (fe) .  P  {K\-:Z^  (fc)  ■  } 

'  ^  p  (^j  (fc)  (fc) ,  ^  {A;r/.:n7 (fe) » } 

nj=l  rii=lmj=l  mi=l 


Now,  as  earlier,  wherever  we  encounter  the  conditioning  event  we  can  replace  ipj  (k) 
with  y”*-’  (k).  Again,  when  we  have  the  event  we  note  that  this  includes  the 

event  ,  hence  the  above  equation  can  be  replaced  by 

^{C.:n7 £>''-'}  = 

p  (vr  (fe)  jk) ,  P  {Kt:n/ (fc) .  } 

1  m  m  lp(^yf^{k)\)^,\t/,^^-Hk),D>‘-^)x' 
Then  we  can  easily  obtain  the  likelihood 

Aj‘{k)  A  p(sr(A:)re.':7,S'^-'(*),o‘-‘) 

=  M{yJ>{k)-,xUk),f”‘{k)) 

=  l2»r7  («:)p’  exp  |-i  (p”'  (fc)  -  xf  (t))'T“'  (fc)-*  (»”'  (k)  -  X?  (*))  } 


where 


(fc)  =  Pl^  (k)  -  p^p  (k)  -  p"^^^  (k)  +  pp  (fe) 


In  the  above 


Pt{k)  A  £{if(fc)x?(fc)'} 

/f'W  &  (fe)'} 

p‘'"'  (fc)  4  E[xUk)s’;'‘(ky] 

p^l’^k)  4  b|p”-1(J:)x?(A:)'} 


(fe)  =  Xi  (fe)  -  x'l  {k) 

(k)  =  Xi  (k)  -  (k) 

and  x^{k)  is  the  most  recent  fused  estimate/prediction  for  target  U  for  the  hypothesis 
being  considered. 
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Hence  combining  the  cases  of  j  =  1,  1  <  j  <  J  and  again  assuming  that  no  two 
resolved  ground  tracks  which  are  due  to  the  same  target  can  be  associated  with  the  same 
propagation  path  gives  the  following  probability  update  for  the  hypothesis 
using  the  track  estimate  of  update  k,  ipj  (fc): 


p  I  \mi..mj 


1  Ml  Mj 

E  E  E  ••  E 

ni=l  nj  =  l  mi=l  mj  =  l 

rhj^Sj 


(15) 


where,  for  j  =  1: 


and  for  1  <  j  <  J: 


Kf{k)  =  Ar’W=p(!/rW|AS‘;»T.-D'‘‘) 
Kf(k)  =  Af-m  &?(!/?■ 

A”'®  A  p(9f  ((!)|A;r':„7, 


The  likelihood  for  both  cases  is  then 


where 


and 


A?'  (fc)  = 


27rT"' 


(fc) 


exp 


-^upikYTpiky^uj^^ik) 


u^yk)  ^  [y^yk)-xUk)) 

Tp{k)  =  Pl^  {k)  -  pI'P  {k)  -  {k)  +  pp  {k) 

ppk)  ^  p{4(fc)xf(A:)'} 

Pp{k)  =  E[yp{k)yp{ky] 

Ppyk)  =  P{x'*(/c)y7^(fc)'} 

Pppk)  =  E[yp{k)x1{ky] 

if  {k)  —  Xi  (k)  -  if  (k) 
yp  (k)  =  Xi  {k)  -  yp  {k) 


(16) 


(17) 


In  the  above,  xf  (k)  is  the  most  recent  fused  estimate/prediction  for  the  target  ti  that 
yj  ^  (k)  is  assumed  to  be  associated  with  for  the  hypothesis  being  considered.  Note  that 
xpk) ,  k  >  0  are  predictions  based  on  fused  data  available  up  to  time  k  —  1,  whereas 
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(k) ,  0,  are  fused  estimates.  The  predictions  {k)  are  used  when  no  fused  estimates 

are  available  for  target  ti  at  time  k. 

Consider  now  the  estimate  x’l^^  (fc)  of  Xi  {k)  and  its  covariance  (k)  obtained  by 
the  fusion  of  x^{k),h  =  0,1,  ..,H  —  1  with  (fc) .  This  is  again  derived  using  the 
fundamental  equations  of  linear  estimation  and  essentially  the  same  reasoning  as  shown 
earlier  for  A:  —  0  . 

The  resulting  fusion  equations  are  hence 

4+‘  (k)  =  i?  (A)  +  (k)  -  (k))  T^‘  (t)-‘  (yj>  (k)  -  if  (k))  (18) 


(k)  =  Pt  (k)  -  [Pt  (k)  -  '  (k))  T"‘  (k)-'  (Pt  (k)  -  P"'"  (k))  (19) 

Some  comments  are  warranted  regarding  a  further  point  of  divergence  from  the  work 
described  in  [17],  [18]  and  [19].  In  those  papers,  a  description  is  given  of  the  application 
of  Gaussian  mixtmre  equations  [4]  for  describing  the  probability  density  function  of  fused 
tracks.  The  idea  is,  given  a  set  of  path  dependent  hypotheses  which  only  differ  in  their  path 
assignments  (ie,  they  all  correspond  to  the  same  path  independent  hypothesis),  one  can 
construct  a  fused  estimate  for  the  corresponding  path  independent  hypothesis  which  is  a 
Gaussian  mbctme  of  the  estimates  for  the  path  dependent  hypotheses.  This  idea  developed 
largely  from  a  limitation  of  the  static  algorithm.  Because  the  hypothesis  probabilities 
in  the  static  algorithm  are  based  on  a  single  update  of  the  tracks,  there  is  usually  not 
enough  information  available  to  achieve  a  single  dominant  path  dependent  hypothesis. 
Usually,  however,  if  one  groups  the  path  dependent  hypotheses  into  their  corresponding 
path  independent  hypotheses,  the  result  is  a  smaller  number  of  hypotheses  with  significant 
probabilities.  This  results  in  greater  ease  of  presentation  and  interpretation  of  data.  This 
approach  does  have  a  down  side,  however.  Knowing  the  correct  path  to  be  associated 
with  a  particular  track  in  radar  coordinates  is  critical  to  accurately  mapping  it  to  ground 
coordinates,  and  hence  obtaining  the  best  fused  estimate.  The  Gaussian  mixture  approach 
as  applied  to  this  problem  does  not  achieve  this.  Instead  it  simply  gives  a  “best”  estimate 
in  a  statistical  sense  based  on  the  relative  probabilities  of  the  propagation  paths.  In  a 
sense,  this  estimate  is  always  wrong,  in  that  it  does  not  correspond  to  any  particular 
path  transformation.  Furthermore  the  relative  probabilities  of  grouped  path  dependent 
hypotheses  in  the  static  algorithm  generally  vary  significantly  from  one  update  to  the  next. 
This  leads  to  undesirable  “jitter”  in  the  fused  tracks  corresponding  to  the  path  independent 
hypotheses.  With  the  introduction  of  the  dynamic  algorithm  this  conundrum  is  largely 
resolved.  Because  the  DMPTF  algorithm  determines  its  hypothesis  probabilities  over 
time,  much  more  data  is  available  for  determining  relative  hypothesis  probabilities.  The 
outcome  is  that  after  the  initiation  of  a  hypothesis  tree,  a  small  number  (often  one)  of  path 
dependent  hypotheses  quickly  become  dominant.  The  result  is  that  there  is  no  need  for 
calculating  Gaussian  mixtures  of  the  path  dependent  hypothesis  estimates  to  determine 
equivalent  path  independent  estimates.  As  a  result  the  aforementioned  approach  was 
discarded. 


29 


DSTO-RR-0223 


6.2  Dealing  with  Changes  in  the  Number  of  Tracks 

A  common  phenomenon  that  occiu'S  is  the  dropping  out  and  initiation  of  multipath 
tracks  as  ionospheric  conditions  change  with  time.  Without  being  able  to  deal  with  this 
phenomenon,  the  DMPTF  algorithm  would  have  little  application.  This  section  deals 
with  additions  to  the  DMPTF  algorithm  which  enable  it  to  cope  with  both  increases  and 
decreases  in  number  of  ground  tracks  with  time. 

First  let  us  consider  what  happens  when  the  number  of  tracks  increases  from  update 
k  to  k  +  1.  The  approach  to  teike  in  this  situation  is  to  predict  the  hypothesis  estimates 
from  update  k  to  update  A:  +  1  exactly  as  for  the  case  when  the  number  of  tracks  does 
not  change,  using  equations  15,  16,  18,  19,  and  then  to  extend  the  tree  using  the  newly 
started  tracks  ,  j  =  1, [J{k  +  1)  -  J  {k)]  .  The  tree  is  extended  in  the  same  way 

as  was  done  during  the  initial  creation  of  the  original  tree,  using  equations  2,  3,  5,  7,  and 

8. 

Now  consider  what  happens  when  the  number  of  tracks  decreases  from  update  k  to 
k  +  1.  Let  us  start  with  a  simple  example  to  help  us  visualize  what  needs  to  be  done 
when  the  number  of  tracks  decreases.  Consider  the  case  of  3  tracks,  2  paths  per  track 
(only  OTHR  tracks)  at  time  k  .  Now  look  at  the  sample  space  fl,  of  all  possible  outcomes 
(where  each  path  dependent  hypothesis  is  a  single  outcome)  as  shown  below.  Note  that 
path  dependent  hypotheses  that  have  two  tracks  with  the  same  propagation  path  for  the 
same  target  are  included  (they  are  not  in  the  actual  implementation). 


^  =  {Allli  -^112)  Ai21,  Ai22,  A123} 


where 


Am  = 

-  f 

- 

\121 

^111) 

\211 

'^111) 

\221 

''ill) 

\112 

''111) 

\122 

''ill) 

\212 

''ill) 

\222 

''ill 

A112  = 

-  f 

-  mii2) 

\  121 
''ll2> 

x211 

Aii2t 

\221 

''112) 

\112 

''112) 

xl22 

''112) 

x212 

''112) 

\222 

''112 

A121  = 

-  f 

-  mi2i> 

\  121 
■^121) 

\211 

^1211 

\221 
''121 ) 

\112 
''121 ) 

\122 

''12I) 

\212 

''12I) 

\222 

''121 

A122  = 

-  \^122) 

^121 

''1221 

x211 

''122) 

x221 

''122) 

>112 

''122) 

Vl22 

''122) 

\212 

''122) 

\222 

''122 

A 123  = 

-  \''123) 

\  121 
^123) 

\211 

''123) 

\221 

''123) 

\112 

''123) 

xl22 

''123) 

\212 

''123) 

\222 

''123 

Now,  the  above  is  only  one  possible  grouping  of  path  dependent  hypotheses.  Let  us 
now  group  all  the  path  dependent  hypotheses  that  have  identical  first  and  last  target 
numbers,  ie,  ignoring  the  target  number  of  the  second  track.  This  is  performed  to  see 
what  needs  to  be  done  if,  say,  track  2  no  longer  exists  at  time  fc  +  1 .  Note  that  we  will 
treat  target  number  3  as  being  the  same  as  target  number  2  for  the  purpose  of  grouping. 
The  reason  that  this  can  be  done  is  that  during  hypothesis  tree  generation  the  target 
numbers  are  incremented  simply  to  signify  that  a  different  target  is  introduced.  Once  the 
target  number  of  the  second  track  becomes  irrelevant,  as  is  the  case  here,  the  third  track 
can  only  be  due  to  the  same  target  or  a  different  target  to  the  first,  and  thus  can  simply 
be  given  a  label  of  1  or  2,  the  former  label  applying  when  it  is  the  same  target  and  the 
latter  when  it  is  not.  We  then  obtain  the  following 

0  =  Ai,(2,3)} 
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The  two  outcomes  listed  on  the  right  hand  side  of  the  above  equations  can  be  decomposed 
into  the  following 


where 


and 


where 


\  _f\l*l  \2*1  \1»2  \2*2\ 
^1*1  — 


\121  \111  \1211 
''1*1  —  I''llli''lll)''l21)''l21/ 
\2*l_r\211  x221  \211  \221\ 
''1*1  —  \''lll>''lll>''l21>''l21/ 

\1*2_  /\112  \122  \112  \122\ 
''1*1  —  t''lll)''lll)''l21)''l21/ 

\2*2_f\212  \222  x212  x222\ 

''1*1  —  I''lll)''lll»''l21»''l21/ 


\  _  / \1*1  \2*1  \1*2  x2*2  \ 

''l*(2,3)  -  \''l*(2,3)>''l*(2,3)>''l*(2,3)'''l*(2,3)/ 


''l*(2,3) 

\2*1 

''l*(2,3) 

\1*2 
^1*(2,3) 


_f\lll  \121  \111  \121  \111  \12n 

—  \''ll2)  ''ll2i  ''122)  ''122)  ''123)  ''l23 / 

f\211  \221  \211  \221  \211  \221  \ 

\-'ll2)  ''112)  ''122)  ''122)  ''123)  ''l23/ 

_  rill2  xl22  \112  \122  \112  \122\ 

—  \  ''112 )  ''112 )  -'122 )  ''122 )  ''123 )  ''l23  / 

\2*2  _r\212  \  222  \212  \222  \212  \  222\ 

''l*{2,3)  “  ■i''ll2)''ll2)''l22)''l22)''l23)''l23J 

Prom  the  above,  it  is  clear  that 

p  {Ai*i}  =  p  {Ai:i } + p  {A?:i } + p  {Ai:f } + p  {a?:?} 


where 


and 


p  { Ai:l }  =  p  { A}il }  +  P  {  A}fl }  +  P  { A“1 }  +  P  { } 

{ A?:l }  -  P  { A?il }  +  P  {  Affl }  +  P  {Xll\ }  +  P  { } 
p  {Ai:?}  =  p  {Ail?}  +  P  {Ai??}  +  p  {Al^?}  +  P  {All?} 

P  {A?:?}  =  P  {A?i?}  +  P  {A???}  +  P  {A?l?}  +  P  {A?l?} 

p  {Ai*(2,3)} = p  {ai:[2,3)} + p  {a?:;2,3)} + p  {a{:?2,3)} + p  {a?:?2.3)} 


where 

P  {a11!(2,3)}  =  P  {Ali^}  +  P  {Al?!}  +  P  {AI22}  +  P  {AI22}  +  P  {Am}  +  P  {Alls} 
p  {^?:(2,3)}  =  P  +  P  +  -p  {^iM}  +  P  {^m}  +  P  {^m}  +  P  {■'la} 

P  {^l*(2,3)}  “  P  {^12}  +  P  {-^112}  +  P  {^122}  +  P  {^122}  +  P  {-^123}  +  P  {^123} 

p  {>f:(2,3)}  =  P  +  P  {Af51}  +  P  {ASl}  +  p  {Afi}  +  P  {Afg}  +  P  {Afi} 
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If  we  now  rename  the  hypotheses  as  follows 


\mi»m3  _  \m1m3 

^ni*n3  min(n3) 

we  have  the  leaves  of  the  hypothesis  tree  for  two  tracks  and  their  corresponding  probabili¬ 
ties.  This  regrouping  can  be  done  at  update  k  once  it  is  known  that  track  2  will  terminate 
at  update  k  +  1  . 

Now  consider  the  target  state  estimates  by  looking  at  some  examples.  Consider  the 
new  hypothesis 


ill  A  \i*i  _  / \111  \121  \111  \12n 
''ll  —  ''l*!  “  I ''ll  lie'll!  >''121  >''121/ 


In  the  above  grouped  hypothesis,  Ajj}  and  Aj^l  have  one  target  estimate  xi  {k)  associated 
with  them,  whereas  Aj^i  and  A}^!  have  2  associated  target  estimates  xi  (k),  X2{k).  The 
new  grouped  hypothesis  should  only  have  target  estimates  that  involve  tracks  1  and  3,  since 
target  estimates  that  only  involve  track  2  must  be  terminated  when  track  2  terminates. 
Now,  in  all  of  Aii},Aif},Aj[2i>'^i2i!  estimate  xi  {k)  involves  tracks  1  and  3  (and  also  track 
2  for  Alii,  A|ii)  hence  xi  (k)  is  retained.  Estimate  X2  {k)  in  A121,  AJ21  only  involves  track 
2  and  hence  is  removed.  Now  we  are  left  with  four  estimates  xi{k) ,  one  for  each  of  the 
path  dependent  hypotheses.  They  are  highly  correlated,  with  the  extent  of  the  correlation 
again  being  extremely  difficult  to  ascertain.  It  is  proposed  that  the  best  approach  here  is  to 
compute  the  equivalent  Gaussian  to  this  Gaussian  mixture  to  form  the  new  xi  (fc).  Their 
high  correlation  is  expected  to  generally  result  in  their  means  being  fairly  close  in  value 
and  hence  the  approximation  in  forming  the  equivalent  Gaussian  should  be  reasonably 
good. 

Now  consider  the  new  hypothesis  AJ2  —  A|*|2  3) 


\11  A  \1*1 
''12  —  ''l*(2,3) 


/illl  il21  3111  \121  %111  \12n 

\''ll2)  ''112)  ''122)  ''122)  ''123)  ''123/ 


In  the  above  grouped  hypothesis 


•  There  are  no  hypotheses  that  have  only  one  target  estimate  associated  with  them. 

•  Hypotheses  A}  12)  A112  A122)  ^122  target  estimates  Xi  (k),  X2  (k)  associated 

with  them.  Neither  of  the  estimates  only  involve  track  2  in  any  of  the  hypotheses, 
hence  both  estimates  are  retained. 

•  Hypotheses  Aj^j  and  Ajls  have  3  associated  target  estimates  x\  {k),  X2{k),  X3{k). 
Estimate  X2  {k)  only  involves  track  2  and  hence  must  be  removed. 


Now  we  are  left  with  six  estimates  xi  (fc);  one  for  each  of  the  path  dependent  hypothe¬ 
ses.  We  then  compute  the  equivalent  Gaussian  to  this  Gaussian  mixture  to  form  the  new 
xi  {k).  We  are  also  left  with  four  estimates  X2  {k) ,  one  for  each  of  A} (2  )  A|i2  )  A122  )  -^122  ) 
and  two  estimates  X3  {k)  ,  one  for  each  of  AI23  and  AI23  ■  Now  estimates  X3  (k)  for  each  of 
A}23  and  A}23  are  associated  with  track  3,  and  estimates  X2  {k)  for  each  of  A}}^  ,  Ajf^  )  A122  ) 
A122  are  also  (at  least  in  part)  associated  with  track  3.  We  hence  compute  the  equivalent 
Gaussian  to  this  Gaussian  mixture  (of  6  Gaussians)  to  form  the  new  X2  {k). 
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We  can  now  generalize  from  the  above  specific  case  to  propose  a  simple  general  algo¬ 
rithm  for  performing  the  regrouping  for  any  number  of  tracks  and  propagation  paths  as 
follows.  Starting  with  a  set  of  path  dependent  hypotheses  where  track  tk 

terminates,  perform  the  following  steps: 

1.  Rename  each  of  {^)  by  removing  subscript  uk  and  superscript  rriK- 

Delete  all  target  estimates  that  only  involve  a  contribution  firom  tk.  Rename  tracks 

TK+1,-,TJ  to  tk,..,tj-i. 

2.  For  each  of  the  newly  named  (k)  (where  subscript  uk  and  superscript 

niK  are  now  removed,  and  tracks  tk+i,..,tj  are  renamed  tk,  rename  the 

hypothesis  using  the  most  compact  notation  that  allows  description  of  tracks  intro¬ 
ducing  new  targets  or  tracks  being  associated  with  previously  hypothesised  targets. 
Reordering  of  subscript  numbers  must  be  performed  for  some  hypotheses  so  that 
the  first  new  target  identifier  is  1,  the  second  is  2,  etc.  Rename  the  target  estimates 
associated  with  each  hypothesis  as  the  corresponding  hypothesis  is  renamed. 

3.  Group  all  hypotheses  that  now  have  identical  labels  (subscripts  and  superscripts)  to 
form  new  merged  hypotheses.  The  probabilities  of  the  new  hypotheses  are  set  to  the 
sum  of  the  probabilities  of  the  hypotheses  from  which  they  are  formed. 

4.  For  each  new  merged  hypothesis,  form  the  equivalent  Gaussian  from  the  (unnor¬ 
malised)  Gaussian  mixture  that  corresponds  to  the  hypotheses  from  which  the  merged 
hypothesis  is  formed. 

If  the  number  of  tracks  decreases  by  more  than  one,  the  above  grouping  is  repeated 
for  each  track  that  has  terminated.  The  hypothesis  tree  is  then  ready  for  the  prediction 
step  from  update  fc  to  fe  -h  1  which  is  performed  as  described  previously. 

The  equations  for  the  replacement  of  the  unnormalised  Gaussian  mixture  associated 
with  a  group  of  hypotheses  by  a  single  Gaussian  probability  density  function  are  given 
by  the  following.  Consider  N  estimates,  Xj,  with  covariances  Pj  and  probabilities  pj,  j  — 
where  the  hypotheses  that  the  estimates  are  associated  with  are  mutually  exclusive 
but  not  exhaustive,  that  is  p  =  Ylf-iPj  ^  1-  The  mean  that  is  used  for  approximating 
the  Gaussian  distribution  is; 


1  ^ 

^  =  (20) 

and  its  associated  covariance  is: 

(21) 

where 

N 

(^j  ~  (®j  “  ^ypj 

j=i 


33 


DSTO-RR-0223 


Note  that  the  factor  1/p  which  appears  in  front  of  the  sum  for  both  the  estimate  and 
covariance  fusion  normalises  the  hypothesis  probabilities  with  respect  to  the  probability 
of  the  merged  hypotheses.  Equations  20  and  21  are  referenced  from  [4]  page  47,  where 
there  is  some  explanation  of  their  derivation  as  well.  Note  that  they  have  been  modified 
to  account  for  the  fact  that  the  hypotheses  are  not  exhaustive. 

To  demonstrate  how  the  above  algorithm  works  in  a  slightly  more  complicated  case, 
consider  the  application  of  the  algorithm  to  two  specific  example  hypotheses  from  a  hy¬ 
pothesis  tree  for  six  radar-coordinate  tracks.  Consider  the  hypotheses  A123123 
with  probabilities  Pa,  and  p;,  respectively.  The  former  has  three  target  estimates,  xi,X2, X3 
associated  with  it,  whereas  the  latter  hypothesis  has  four,  ie  ii,X2,X3,X4  .  Let  track  2 
terminate.  After  step  1  the  new  hypothesis  names  are  Ajm®  and  AJ3143  respectively.  For 
the  former  all  three  target  estimates  xi,X2,X3  remain,  whereas  for  the  later  the  estimate 
X2  is  removed  leaving  xi,  X3,  X4.  After  step  2  the  two  hypotheses  have  identical  hypothesis 
names  and  estimate  labels,  ie  Ajim  and  xi,X2,X3  respectively.  Hence  the  two  hypotheses 
are  then  merged  to  form  the  new  hypothesis  AJ2132  with  probability  p  =  Pa+  Pb  ■  The 
formation  of  the  three  Gaussian  mixtures  and  their  equivalent  Gaussians  is  performed  by 
a  straight  forward  application  of  equations  20  and  21. 


6.3  Track  Dependence 

The  likelihood  and  fused  estimate  equations  in  both  the  static  and  dynamic  cases 
incorporate  cross-covariance  terms  so  that,  in  theory,  dependence  between  estimates  rpjik) 
can  be  accommodated.  The  main  causes  of  dependence  between  tracks  are 

•  common  process  noise  [5], 

•  tracks  being  updated  using  common  measurements  [19],  [7],  [13],  and 

•  common  propagation  path  segments  [19]. 

In  practice,  however,  it  was  found  that  the  cross-covariance  terms  are  extremely  dif¬ 
ficult  to  determine  because  of  the  complexity  of  the  interactions  involved.  Research  was 
undertaken  to  determine  the  cross-covariance  terms,  or  at  lesist  to  approximate  them,  with 
at  this  stage  little  success. 

The  approaches  for  dealing  with  dependence  between  tracks  that  were  looked  at  were 

1.  Assume  independence,  ie,  set  cross-covariances  to  zero. 

2.  Directly  calculate  cross  covariances. 

3.  Simple  approximations  to  cross-covariances,  based  on  judgement. 

4.  The  covariance  intersection  algorithm  [27]. 

5.  Rearrangement  of  equations  to  reconstruct  measurements. 

6.  Using  the  measurements  from  which  the  original  unfused  tracks  have  been  formed. 
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A  summary  of  the  approaches  and  the  results  obtained  upon  testing  is  given  in  the  following 
paragraphs. 

The  first  approach,  which  is  also  the  easiest,  is  to  simply  ignore  the  interdependencies 
between  the  multipath  tracks,  ie,  to  set  the  cross-covariance  terms  in  the  likelihood  cal¬ 
culations  and  fusion  equations  to  zero.  This  was  tested  and  used  as  a  reference  against 
which  the  performance  of  the  other  approaches  was  compared.  Somewhat  surprisingly, 
while  obviously  not  optimal,  the  approach  was  found  to  provide  reasonable  fusion  and 
likelihood  calculations,  with  the  major  negative  effect  being  that  the  covariances  of  the 
fused  estimates  shrink  somewhat  too  quickly  over  time.  This  was  found  to  result  in  a 
particular  favoured  hypothesis  having  its  probability  approach  1  more  rapidly  than  would 
be  expected.  It  also  led  to  the  resulting  estimates  being  very  “smooth” ,  when  compared 
with  the  original  tracks,  again  more  so  than  an  optimal  algorithm  would  be  expected  to 
produce. 

The  second  approach  as  suggested  in  [19],  and  described  in  [5]  is  extremely  unwieldy 
when  applied  to  DMPTF,  because  of  the  number  of  tracks  that  need  to  be  fused  and  the 
asynchronicity  that  must  be  accommodated.  This  technique  also  puts  restrictions  on  the 
models  used  by  the  trackers  that  feed  information  to  the  DMPTF  algorithm.  Additionally, 
the  approach  would,  if  in  fact  it  could  be  implemented,  at  best,  only  solve  the  process  noise 
aspect  of  the  dependence.  For  these  reasons,  this  technique  was  not  considered  further. 

The  third  approach  is  to  provide  a  simple  approximation  for  the  cross-covariance  ma¬ 
trices,  as  suggested  in  [19],  and  [5].  Considerable  experimentation  was  performed  with 
this  technique,  but  it  was  found  that  this  approach  simply  did  not  work  for  DMPTF.  A 
geometric  combination  [5]  p455  of  the  covariances  of  the  two  dependent  estimates  as  well 
as  several  other  types  of  approximation  of  the  cross-covariance  terms  were  tried,  none  of 
which  came  even  close  to  giving  satisfactory  performance.  The  resulting  covariance  sums 
(shown  in  equations  4  and  17)  would  sooner  or  later  result  in  a  non-singular  covariance 
matrix,  leading  the  DMPTF  software  execution  to  terminate.  The  only  way  that  this 
could  be  avoided  was  by  setting  the  cross-covariance  terms  to  values  very  close  to  zero, 
effectively  approximating  the  first  approach  listed  above.  The  suspected  reason  is  that  the 
dependence  between  new  data  and  previously  fused  data  varied  substantially  with  time 
and  between  hypotheses,  making  it  impossible  for  a  simple  cross-covariance  approximation 
to  be  used  successfully. 

The  fourth  option  considered  is  the  covariance  intersection  (Cl)  algorithm  [27]  (also 
referred  to  as  Gaussian  Intersection).  The  aim  of  this  approach  is  to  avoid  the  need  to 
calculate  or  estimate  dependence  between  data  which  is  to  be  fused,  while  stiU  arriving 
at  a  result  that  is  optimal  in  some  sense.  A  short  summary  of  the  algorithm  and  some 
experimental  findings  will  be  given  here;  full  details  of  the  technique  can  be  found  in  [27]. 

Consider  two  estimates,  Xa  and  Xb  with  means  a  and  b  and  covariances  Paa  and  Pbb 
respectively.  The  covariance  intersection  algorithm  forms  the  fused  estimate  c  with  co- 
variance  Pec  by  using  the  following  equations 

+  (22) 

where  OJ  €  [0, 1] 

The  free  parameter  lo  can  be  used  to  optimize  the  fused  estimate  with  respect  to  different 
performance  criteria  such  as  minimizing  the  trace  or  the  determinant  of  Pcc.  The  result. 
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Figure  6:  Example  of  fusion  using  the  Covariance  Intersection  Technique 


it  is  argued,  is  a  conservative  and  non-divergent  update  of  the  fused  estimate.  Figure 
6  shows  a  simple  example  of  the  fused  estimate  produced  by  the  covariance  intersection 
technique.  Three-sigma  ellipses  for  two  coincident  estimates  are  shown  (labelled  Pa  and 
P(,),  as  well  as  two  ellipses  representing  the  fused  estimate.  The  Ellipse  labelled  “Indep.” 
represents  the  fused  estimate  obtained  by  assuming  that  the  estimates  Xa  and  are 
independent.  The  ellipse  labelled  “Cl”  represents  the  fused  estimate  obtained  using  the 
covariance  intersection  technique  with  determinant  minimization. 

The  Cl  technique  looked  quite  promising  initially;  however,  in  testing  done  to  date  its 
performance  has  been  disappointing.  Updating  of  estimates  with  u>  set  to  minimize  the 
trace  as  well  as  the  determinant  of  the  fused  estimate’s  covariance  (within  the  constraints 
of  equations  22)  was  implemented  in  the  DMPTF  algorithm,  and  tested.  The  outcome  was 
that  when  determinant  minimisation  was  used,  the  Cl  algorithm  would  usually  either  select 
one  or  other  of  the  two  estimates  (ie,  a;  w  0,  or  a;  «  1)  as  its  fused  estimate,  and  tended  to 
jump  between  the  two  in  subsequent  updates.  Using  trace  minimization  provided  a  more 
stable  result,  but  once  again  the  algorithm  would  choose  one  of  the  associated  estimates  as 
the  fused  estimate,  rather  than  providing  a  “weighted  average”  of  the  two.  While  testing 
performed  to  date  has  shown  poor  results,  leading  to  the  technique  being  put  aside,  further 
investigation  to  analyse  the  reasons  for  the  unsatisfactory  performance  may  be  worthwhile, 
and  possibly  lead  to  improvements  in  performance. 
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The  fifth  option  considered  was  to  perform  measurement  reconstruction  [15],  [16],  by 
manipulation  of  the  Kalman  filter  equations  which  are  used  to  perform  tracking.  This 
approach  offers  the  opportunity  to  account  for  the  effects  of  common  process  noise  with¬ 
out  calculating  the  cross-covariances.  It  would  not  account  for  the  effects  of  tracks  being 
updated  with  common  measurements  or  the  common  propagation  path-segment  effects. 
For  our  application,  it  was  found  to  offer  no  advantages  over  simply  using  the  measure¬ 
ments  themselves,  since  using  the  measurements  would  be  no  more  difficult  than  accessing 
the  information  required  for  measurement  reconstruction.  Measurement  reconstruction 
also  has  the  disadvantage  of  not  being  able  to  account  for  common  measurements  as  can 
obviously  be  done  if  the  actual  measurements  are  used  in  the  fusion  process,  and  is  more 
difficult  to  generalize  to  a  multi-sensor  application  where  each  sensor  could  potentially  use 
a  different  Kalman  filter  (or  some  other  filter  for  that  matter).  As  a  result  this  approach 
was  not  investigated  further. 

The  sixth  approach  that  was  considered  was  to  associate  and  fuse  the  measurements 
from  which  the  radar  coordinate  tracks  were  formed.  It  should  be  emphasised  here  that 
only  the  measurements  which  have  already  been  selected  by  the  tracking  algorithm/s  are 
considered,  thus  avoiding  the  computational  load  associated  with  the  centralized  tracking 
approach.  Doing  this  immediately  avoids  the  dependence  due  to  common  process  noise, 
and  offers  the  potential  to  remove  the  effects  of  common  measurements  by  simply  inspect¬ 
ing  the  measurements  prior  to  fusion.  This  approach  still  does  not  address  the  common 
propagation  path  segment  contribution  to  track  dependence,  notwithstanding  that  this 
contribution  may  be  less  difficult  to  deal  with  once  the  other  effects  are  removed.  The 
approach  has  considerable  merit;  hence  the  association  and  fusion  equations  required  to 
use  measurements  have  been  derived  and  are  presented  later  in  this  report.  The  approach 
has  not  as  yet  been  implemented  for  testing  however. 

To  date,  the  first  option,  ie,  assume  that  the  data  is  independent,  is  what  is  imple¬ 
mented  in  the  DMPTF  algorithm. 
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7  Performance  Testing  of  DMPTF  Algorithm 

The  DMPTF  algorithm  has  been  implemented  in  C++  and  integrated  into  a  devel¬ 
opment  test-bed  in  order  to  determine  its  performance  and  assess  further  development 
requirements. 

For  the  initial  testing  of  the  DMPTF  algorithm,  pruning  was  set  to  retain  the  100 
highest  probability  hypotheses,  and  the  cross-covariance  terms  of  equations  4  and  17  were 
set  to  zero.  Stored  examples  of  real  data  from  the  Jindalee  OTHR  at  Alice  Springs, 
Australia  were  used  as  inputs  to  the  algorithm  and  performance  ascertained.  Two  types 
of  ionospheric  coordinate  registration  (CR)  data  were  used  for  the  path  transformations. 
Actual  CR  data  based  on  ionospheric  sounder  measurements  was  used,  as  well  as  data 
based  on  a  simplified  two-layer  ionospheric  model.  The  reason  for  using  the  simplified 
model  in  some  of  the  testing  was  to  enable  control  over  more  parameters. 

During  the  early  testing  it  became  evident  that  performance  could  be  improved  by 
making  some  minor  modifications  to  the  DMPTF  algorithm.  It  should  be  noted  that  the 
changes  to  the  algorithm  were  essentially  of  an  empirical  nature,  giving  improved  perfor¬ 
mance  for  the  application  in  question,  but  not  adding  to  the  theory.  The  first  modification 
was  made  because  in  some  cases,  if  a  hypothesis  which  was  originally  of  highest  probability 
became  inappropriate  several  updates  later,  the  algorithm  would  sometimes  not  readjust 
sufficiently  quickly.  The  cause  of  this  was  that  the  probabilities  of  non-favoured  hypothe¬ 
ses  would  become  extremely  small  after  several  updates,  thus  making  it  difficult  for  them 
to  quickly  become  the  highest  probability  hypothesis  when  conditions  changed.  The  rem¬ 
edy  for  this  was  to  set  a  minimum  value  of  the  probability  for  each  hypothesis.  At  every 
update  the  hypothesis  probabilities  would  be  calculated  using  the  equations  outlined  in 
the  earlier  sections,  then  for  all  hypotheses  (among  the  top  100)  that  had  a  probability 
below  the  set  threshold,  their  probabilities  would  be  set  to  the  value  of  the  threshold,  after 
which  the  probabilities  of  all  the  hypotheses  would  be  re-normalized.  The  effect  of  this 
change  was  to  allow  the  algorithm  to  adapt  more  quickly  to  changing  input  data.  Some 
“tuning”  of  this  modification  is  required  to  obtain  best  performance. 

The  second  modification  made  was  to  change  the  way  that  hypothesis  probability  up¬ 
dating  was  performed  for  hypotheses  containing  targets  with  a  single  contributing  track 
(single-track-per-target  hypotheses).  This  modification  was  required  for  reasons  that  were 
essentially  an  outcome  of  limitations  in  the  quality  of  the  ionospheric  and  track  data,  as 
well  as  limitations  in  the  modeling  of  the  system.  The  major  effect  that  caused  problems 
was  that  for  the  single-track-per-target  hypotheses,  where  the  only  difference  is  the  prop¬ 
agation  path  assignment,  after  several  updates  one  hypothesis  would  become  much  more 
probable  than  all  of  the  others.  If  the  data  and  models  were  a  perfect  representation  of 
reality,  this  would  not  be  a  problem.  However,  because  this  is  not  the  case,  essentially  se¬ 
lecting  one  of  the  contending  hypotheses  has  often  not  been  appropriate.  The  adjustment 
made  was  to  use  a  modified  likelihood  update  equation  in  place  of  equation  16.  The  re¬ 
placement  equation  at  this  stage  is  equation  3  and  is  essentially  an  ad  hoc  modification.  In 
practice  this  change  makes  the  single-track-per-target  hypotheses  “fail-back”  hypotheses, 
whose  relative  probability  with  respect  to  one  another  does  not  change  with  time  (but  may 
change  with  respect  to  other  hypotheses) ,  and  whose  probability  will  only  become  close  to 
1  if  all  other  hypotheses  are  a  poor  match  to  the  data.  A  more  appropriate  modification  to 
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equation  16  may  be  found  with  further  study,  although  the  present  modification  appears 
adequate  from  a  performance  point  of  view. 

As  an  indication  of  algorithm  performance,  Figures  7  and  8  show  a  typical  multipath 
track  scenario  and  the  results  of  fusion  using  DMPTF  over  a  total  of  50  updates.  Figure  7 
shows  track  data  received  from  the  OTHR  in  radar  coordinates  (slant  range  and  azimuth). 
Figure  8  shows  the  fused  output  of  DMPTF  in  ground  coordinates  (range  and  azimuth) 
using  the  simplified  two-layer  ionospheric  model,  resulting  in  four  possible  propagation 
paths  for  each  track.  In  the  ionospheric  model,  the  height  of  the  E  layer  was  set  to  100 
km  with  a  standard  deviation  of  10  km  and  the  F  layer  height  was  set  to  300  km  with  a 
standard  deviation  of  20  km.  The  underlying  truth  of  this  example  is  represented  by  Figure 
9,  which  shows  track  data  received  from  a  microwave  radar  for  two  targets  travelling  from 
the  bottom  right  to  the  top  left  of  the  figure.  Before  continuing,  let  us  also  define  the  two 
letter  notation  used  in  the  following  text  to  denote  the  highest  probability  propagation 
path  assignment  made  by  the  DMPTF  algorithm.  The  first  letter  denotes  the  ionospheric 
path  taken  by  the  transmitted  signal  (ie,  from  the  radar  to  the  target)  and  the  second 
letter  denotes  the  retmn  path  (ie,  from  the  target  to  the  radar).  For  example  an  EF 
propagation  path  assignment  represents  transmission  via  the  ionospheric  E  layer  and  a 
return  signal  propagating  via  the  F  layer. 

The  input  data  to  the  DMPTF  algorithm  (Figure  7)  consists  of  seven  tracks  in  total. 
Tracks  1,  2  and  3  are  due  to  target  1  in  Figure  9  and  tracks  4,  5,  6  and  7  due  to  target  2. 
Note  that  the  seven  tracks  do  not  all  exist  for  the  entire  duration  of  the  example.  TVacks 
1,  2,  4  and  5  exist  for  the  entire  duration  while  track  3  exists  at  the  beginning  but  drops 
out  towards  the  end.  TFack  7  starts  a  little  way  into  the  example,  and  track  6  starts  a 
little  later  still;  both  tracks  then  run  to  the  end  of  the  duration  of  the  example. 

Figure  8  displays  the  fused  target  estimates  of  the  highest  probability  path-dependent 
hypothesis  at  each  update;  hence  the  figure  does  not  show  the  same  hypothesis  for  all 
updates.  The  first  feature  that  should  be  apparent  about  the  result  is  that  there  are  four 
distinct  tracks  shown  in  the  figure,  describing  the  two  targets  that  actually  exist.  The 
example  illustrates,  as  will  be  explained  below,  that  with  typical  real  data  the  algorithm 
may  not  always  deliver  a  perfect  result,  but  the  results  that  DMPTF  does  give  are  generally 
consistent  with  the  sometimes  limited  and  ambiguous  data  that  is  available.  The  example 
was,  in  fact,  chosen  because  it  highlighted  some  of  the  ambiguities  that  occur  in  common 
scenarios,  which  can  sometimes  lead  to  more  than  one  “correct”  track  to  target  assignment. 
Note  also  that  the  azimuth  scale  has  been  greatly  expanded  relative  to  the  range  scale  to 
highlight  differences  between  target  positions  when  path  assignments  change.  The  scales 
have  not  been  marked  on  the  axes,  however,  to  maintain  an  “unclassified”  classification 
for  this  document. 

Target  1  in  Figure  8  is  the  result  of  fusing  tracks  1,  2  and  3  from  Figure  7,  the 
propagation  path  assignments  being  FF,  EF,  and  EE  respectively.  This  target  is  consistent 
with  the  underlying  truth  for  the  example,  ie,  it  matches  up  with  target  1  of  Figure  9.  This 
consistency  is  maintained  throughout  the  duration  of  the  example,  despite  the  switching  of 
hypotheses  associated  with  changes  in  propagation  path  assignments  for  tracks  associated 
with  target  2,  as  well  as  changes  in  the  total  number  of  tracks  in  the  hypothesis  tree. 

Target  2  in  Figure  8  starts  in  the  bottom  right-hand  corner  of  the  figure  at  the  same 
time  as  target  1.  TVacks  4  and  5,  which  are  the  only  tracks  that  exist  at  this  stage  that  are 
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Figure  7:  OTHR  tracks  in  radar  coordinates 
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actually  due  to  target  2,  are  fused  together  at  the  outset.  The  propagation  assignments 
for  tracks  4  and  5  are  initially  FE  and  EE,  respectively.  As  becomes  evident  later  in  the 
example,  these  initial  propagation  path  assignments  are  incorrect.  The  reason  for  this  is 
that  there  are  only  two  tracks  available  for  association  and  fusion,  resulting  in  ambiguity 
regarding  which  pair  of  propagation  path  assignments  is  appropriate  for  the  target.  This 
ambiguity  is  not  due  to  the  DMPTF  algorithm  itself,  but  results  from  the  geometry  of  the 
propagation  paths.  For  example,  the  relative  positions  of  two  tracks  which  are  due  to  a 
single  target,  but  propagating  via  the  EE  path  and  FE  path  respectively  are  very  similar 
to  those  of  two  tracks  propagating  via  EF  and  FF  respectively.  At  the  points  labelled  A 
and  B  a  third  track  is  initiated,  track  7,  which  is  originally  determined  to  be  a  third  target, 
(and  assigned  the  EE  path),  but  should  be  associated  with  tracks  4  and  5,  ie,  target  2. 
This  incorrect  target  assignment  is  due  to  the  fact  that  the  path  assignments  for  target 
2  are  initially  incorrect.  Due  to  the  incorrect  path  assignments,  the  third  track  cannot 
initially  be  associated  with  target  2  in  the  highest  probability  hypothesis,  so  it  as  assigned 
to  a  third  target.  However,  within  two  updates,  the  hypothesis  that  fuses  tracks  4,  5  and 
7,  with  path  assignments  of  FF,  EF  and  EE  respectively,  becomes  the  highest  probability 
hypothesis,  bringing  the  separate  targets  at  A  and  B  to  C.  The  path  assignments  for  both 
targets  1  and  2  are  then  correct.  When  track  6  is  initiated  at  a  later  time,  it  is  seamlessly 
fused  into  target  2  with  path  assignment  FE,  and  when  track  3  terminates,  leaving  target 
1  with  two  tracks  contributing  to  it  and  target  2  with  four  tracks  contributing  to  it,  the 
target  estimates  continue  smoothly. 

The  example  gives  an  indication  of  the  difficulties  associated  with  the  environment  in 
which  the  DMPTF  algorithm  must  operate.  The  ambiguity  between  path  assignments  can 
sometimes  result  in  an  incorrect,  but  plausible,  set  of  paths  being  assigned  to  a  track.  The 
effect  of  implementing  the  lower  limit  of  probability  for  each  hypothesis  is  also  evident 
in  this  example,  in  that  an  originally  incorrect  assignment  is  rapidly  resolved  when  more 
information,  in  the  form  of  a  third  track,  becomes  available.  Where  the  assignments  are 
correct,  the  example  shows  that  the  DMPTF  algorithm  produces  smooth,  fused  estimates 
of  the  targets  in  their  correct  spatial  positions,  despite  losing  and  gaining  track  information 
as  the  example  evolves  over  time. 

In  addition  to  implementation  in  the  development  test-bed,  the  DMPTF  algorithm  has 
b^n  integrated  into  the  JFAS  radar  developmental  code  for  more  thorough  testing.  To 
this  end  the  radar’s  graphics  displays  were  also  modified  to  enable  display  of  fused  tracks 
that  are  produced  by  the  DMPTF  algorithm.  Recently  an  extended  trial  of  the  algorithm 
has  been  performed  on  the  JFAS  radar,  and  a  large  amount  of  data  collected  for  statistical 
analysis.  While  the  statistical  analysis  has  not  been  completed  yet,  early  indications  are 
that  the  algorithm  performs  well  when  good  coordinate  registration  (CR)  data  is  available. 
When  the  CR  data  is  poor,  as  might  be  expected,  incorrect  associations  can  occur.  Based 
on  the  testing  to  date,  indications  are  that,  at  the  present  level  of  development  of  the 
DMPTF  algorithm  and  the  CR  system,  the  DMPTF  algorithm  is  capable  of  being  used 
in  a  semi-automated  fashion  to  give  advice  to  an  operator,  but  operator  overseeing  is  still 
required.  However,  higher  levels  of  automation  can  be  expected  with  further  development. 
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8  Track  Association  and  Fusion  using 

Measurements 


There  are  important  advantages  in  associating  and  fusing  the  measurements  from  which 
the  tracker  forms  its  tracks  instead  of  fusing  the  tracker’s  estimates,  ie,  to  use  the  tracker 
as  a  “selector”  of  measurements  for  the  DMPTF  algorithm.  This  approach  provides  key 
advantages  of  using  measurements  while  still  relieving  DMPTF  of  the  computational  com¬ 
plexity  associated  with  using  all  available  measurements  in  the  association  and  fusion 
process  (including  the  responsibility  of  track  initiation)  as  was  done  in  [22].  The  ad¬ 
vantages  come  about  from  the  reduced  dependence  of  data  that  is  used  by  the  DMPTF 
algorithm  when  measurements  are  used.  As  was  mentioned  earlier  in  the  report,  track 
dependence  occurs  as  a  result  of  the  following: 


1.  common  process  noise, 

2.  tracks  being  updated  using  common  measurements,  and 

3.  common  propagation  path  segments. 


Of  the  three  causes  of  dependence  above,  the  first  does  not  occur  at  all  when  mea¬ 
surements  are  used  in  place  of  track  estimates.  The  second  source  of  dependence  can 
reasonably  readily  be  removed  by  checking  for  common  measurements  during  the  fusion 
process  and  discarding  multiple  occurences  of  measurements  prior  to  fusion.  The  third 
source  of  dependence  still  remains,  and  must  be  dealt  with  by  other  means.  Hence,  by  fus¬ 
ing  the  measurements  from  which  the  tracker  forms  its  tracks,  two  of  the  three  sources  of 
data  dependence  can  be  removed  with  no  significant  increase  in  computational  complexity. 

Let  us  now  derive  the  association  and  fusion  equations  for  the  case  where  measurements 
are  used.  Note  that  it  is  quite  similar  to  the  derivation  for  the  case  where  estimates  are 
used  which  was  given  in  earlier  sections;  the  differences  come  in  the  details  of  the  equations. 

Consider  the  creation  of  the  hypothesis  tree  and  calculation  of  the  hypothesis  proba¬ 
bilities  and  fused  estimates.  As  indicated  above,  the  measurements  that  are  used  are  those 
that  the  tracking  algorithm  uses  to  produce  the  estimates  ipj  (0) ,  ipj  (1) ,  '4>j  (2) , ...  For  the 
purpose  of  the  subsequent  discussion,  it  will  be  assumed  the  tracking  algorithm  is  a  nearest- 
neighbour  (NN)  tracker  [3].  The  approach  described  below  can  easily  be  generalized  for 
use  with  a  PDA  tracker.  Now  let  us  denote  the  measurements  as  5j  (0)  ,6j  (1)  ,5j  (2)  ,.. 
Note  that  because  of  the  assumption  of  a  NN  tracker  there  is  only  one  measurement  per 
update,  ie,  for  each  estimate  in  slant  coordinates  there  is  only  one  associated  measurement. 
Note  firstly  the  difference  between  the  track  estimate  vector  and  the  measurement  vector. 
The  track  estimate  vector  in  slant  coordinates  is 


Rj 

Rj 

Aj 

Ai 
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The  measurement  vector  in  slant  coordinates  is: 


Consider  first  the  creation  of  the  hypothesis  tree  at  time  k  =  0.  The  derivations  will 
follow  a  very  similar  line  to  that  for  the  case  where  track  estimates  are  used.  Note  also 
that,  as  was  done  in  section  2  the  time  index  (0)  will  be  omitted  for  time  k  =  0.  Prom  the 
first  measurement,  (5i,  the  propagation  path  transformations  can  be  used  to  calculate  the 
corresponding  measurement  in  ground  coordinates  and  its  covariance  for  each 
propagation  path  mi  —  l,..,Mi.  Under  the  assumption  that  there  are  no  false  tracks, 
there  is  only  one  target  association  possible,  namely  target  1.  The  probability  of  each  of 
the  path  dependent  hypotheses  is  then: 

(23) 

where  0^^  =  P  |  and  6^^  denotes  the  event  that  track  Tj  is  propagating  via  propa¬ 
gation  path  mj.  As  for  the  case  of  using  track  estimates,  the  0^’  are  prior  probabilities 
which  are  estimated  using  physical  measurements  of  the  ionosphere. 

For  each  subsequent  measurement,  <5j,  j  =  2, ..,  J,  the  following  recursive  equation  for 
calculating  the  probability  of  from  can  be  derived  in  a  similar  fashion 

to  that  used  when  track  estimates  were  considered: 

nK=l  _mK=l! 
rriK^SK 

where 

KT 

Ar 

and  A'^  =  (Ji, ..,  6k 

If  riK  0  rij  for  all  j  —  1, ..,  K  —  1,  that  is  5k  represents  a  new  target,  there  is  no  prior 
information  regarding  the  state  of  .  Hence  the  likelihood,  ,  is  given  by 

Ar  =  0  (26) 

ym 

where  Vm  is  in  this  case  the  “volume”  in  measurement  space  in  which  the  target  may  be. 

If  nK  —  nj  for  one  or  more  of  j  =  1, ..,  K  —  1,  then  6k  represents  the  same  target  as 
at  least  one  previously  hypothesized  target.  In  this  case  the  new  information  contained 
in  6j  must  be  fused  with  the  previous  estimates.  Assume  that  the  hypothesis 
is  associated  with  T  targets,  1  <  T  <  AT  -  1,  that  the  target  corresponding  to  the  track 
tk  is  ti,  1  <i  <T,  and  that  prior  to  considering  2]^*^  there  have  been  h  track  estimates 


—  P{Zk  1-^1  >••>%-!  i^ni..nK  ) 

A  jriK-i  N 

Py'K  fi^K-l  » ^ni..nif_inK  J 


Argr-P{-C.‘:»"-:‘ 

Mfc 

E  E  Ar/?r 


45 


DSTO-RR-0223 


associated  with  target  ti  for  the  current  time  instant.  Let  the  previous  estimate  of  the  state 
of  target  ti  be  and  its  covariance  Since  it  is  assumed  that  (with  covariance 
is  from  target  ti  the  best  prediction  (the  expected  value)  of  2^''  is 


Eb 


.■niK 

■K  ^1  : 


rriK-i 

‘A'-l 


where 


4  =  Hx^ 


The  covariance  of  —  z^)  is: 


=  E 


^rriK  _  =h 
H 


J^K  _ 


=  E  I  -  Hxi)  -  H  (if  -  Xi)]  [(2^^  -  Hxi)  -  H  (if  - 

Let  if  =  Xj  -  if  and  2^^  =  Hxi  —  2^^  ,  where  Xi  is  the  true  state  of  t«irget  U  ,  then 

fr  . 

E- 


ST 


L/if  -  zT 


H~x^-z^x 


=  £|;f4(x'-)'«'-ffi.''(fr)'-5r  (4)'«'+sr  k")' 
=  £{«x*  (x'“)'h'|  -  (*0’'''} 


Now 


Pt  =  E{x\ 


and  let 


^  E 


{ifK")'} 


ErzT 


Then 


ST  ^  +  jimK 


The  likelihood  of  zT  (assuming  a  Gaussian  distribution)  is  then  given  by 
A^x  ^  |27rS^^|"=exp 


-^^i^T)'{sTr'^T 


(26) 


(27) 


where 


rriK  A  I 


{‘T  -  4) 
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Consider  now  the  estimate  of  Xi  and  its  covariance  obtained  by  the  fusion 
of  x^,h  =  1,2,.., if  —  1  with  .  To  derive  this  we  shall  again  use  the  fundamental 
equations  of  linear  estimation  as  shown  on  pages  44  and  125  of  [4]  and  reproduced  in 
equations  6  in  section  2.  Hence,  replacing  the  terms  in  the  fundamental  equations  of 
linear  estimation  with  their  equivalent  terms  in  this  problem,  and  following  the  procedure 
used  in  section  2  to  derive  equations  7  and  8,  results  in  the  following  fusion  equations: 


=  i?  +  [pf'H'  -  C/f”")  (zS*  -  zfj  (28) 

ph+i  =  _  (pkg'  _  (nff  -  M  (29) 


Now  let  us  consider  updates  k  =  1,2, 3, ...  .  As  before,  the  hypothesis  tree  built  up  at 
k  —  0  can  be  used  as  the  basis  for  the  tree  at  subsequent  updates.  The  derivations  again 
follow  a  similar  line  to  those  for  the  case  of  using  track  estimates.  In  order  to  shorten 
notation  the  symbol  (fc)  is  defined  such  that  {k)  =  (fc) ,  ..,(5j  (k). 

In  the  creation  of  the  hypothesis  tree  (at  update  fc  =  0)  we  derived  the  probability 
P  (0)  }  of  each  path  dependent  hypothesis  nij  £  {1, 2,  ■.,Mj}  ,  rij  € 

{l,2, Bni..nj-i}  given  all  the  data  available  at  update  A:  =  0 ,  ie  the  measurements 
(5i  (0)  ,..,dj  (0)  .  Now  let  us  define  the  following  shorthand  notation 

D'^  -  <5i  (0)  ,..,dj  (0) ,  (5i  (1) , ..,  (1) , ....,  (A:  -  1) , ..,  (A:  -  1) ,  (5i  (A:) , ..,  (A:) 

=  A^(0),A''(1),..,A-'(A:) 


ie,  D*’  represents  all  the  track  data,  in  radar  coordinates,  (measurements  in  this  case) 
available  up  to  and  including  update  k  .  Hence 


P 


pk  }  4  P  I A-'  (0) ,  A-'  (1) , A-'  (k) } 


Consider  the  updating  of  the  probability  of  the  hypothesis  when  the  first 

measurement  for  update  k ,  (ie,  5i  (A:))  is  considered.  Using  Bayes  rule  we  have 


P  {Si  (fe)  P  } 

p{5i{k)\Dk-^) 

p  {^1  jk)  P  } 


Mj  Ml 
nj=l  ni=lmj=l  m\=l 


Now,  as  noted  earlier,  in  the  implementation  of  MPTF  the  ground  co-ordinate  mea¬ 
surements  ^  (k)  are  used  to  perform  the  probability  calculations.  Again,  wherever  we 
encounter  the  conditioning  event  we  replace  Sj  (k)  with  (k).  Now,  when  we  have 
the  event  J^ote  that  this  includes  the  event  .  Hence  the  above  equation 
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can  be  replaced  by 

p{Ki::T/  I'Si  = 

_ p  (k)  p  I-d*~‘  } _ 

\  Mj  M\  _ 

E  E  E  ..  E 

nj=l  ni  =  l  mj  =  l  mi=l 

To  calculate  p  {k)  ,  D^~^)  we  can  again  use  information  regarding  the 

dynamics  of  the  target  that  (k)  is  associated  with.  Now,  as  in  section  6.1,  let  us 
assume  that  the  motion  of  target  U  satisfies  equation  13.  Hence,  given  the  fused  estimate 
xf  {k  -  1)  for  the  hypothesis  being  considered,  where  H  is  the  number  of  track  estimates 
that  were  combined  at  time  —  1  to  obtain  xf  {k  —  1) ,  we  have  (as  in  equations  14) 

{k)  =  F{k-l)  xf  {k  -  1) 

W  =  f  (t  -  1) (t  -  1)  F  (t  -  1)' + 

Using  the  above  we  can  easily  obtain  the  likelihood 

=  |2,rSJ”'  (fc)|-iexp  |-i  (k)  -  (k))' (*)^>  (x”'  (k)  -  i”  {P))| 

where 

{k)  =  H  (k)  Pf  (k)  H  {k)'  -  H  {k)  f/^>  (k)  -  (k)  H  (k)'  +  (k) 

with 

P«{k)  A  F{i»(fe)xP(*)'} 

RT'{k)  =  F{2f’‘(A:)2f‘{A;)'} 

U^‘{k)  A  E{x»(*:)i”'(*:)'} 

[/,”■"(*)  A  E{z"'{k}xUky} 

and 

(A:)  =  Xi  {k)  x^  (k) 
z^^{k)  =  Hx^{k)-z^^  {k) 

Consider  now,  the  updating  of  the  probability  of  the  hypothesis  when  the 

measurement  for  update  k  ,  (ie,  6j  {k) ,  1  <  j  <  J  )  is  considered.  Using  Bayes  rule  we 
have 

p jk)  {k),D^-^) p {XZT/  |a^-^ {k),D’^-^ } 

p{Sj{k)\A^-^{k),D>^-^) 

P  {Sj  jk)  (fc)  p  |A^-1  (fc)  } 

••  Z^  ••  Z^  p  { \  Aj-1  {i.\  nfc-l  1 

nj  =  l  ni  =  lmj  =  l  mi=l  L  f 
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Now,  as  earlier,  wherever  we  encounter  the  conditioning  event  9^^  we  can  replace  5j  (k) 
with  {k).  Again,  when  we  have  the  event  J^ote  that  this  includes  the 

event  9^^  ,  hence  the  above  equation  can  be  replaced  by 

p  (*=)  (k)  ,!)*-■)  p  |aj-‘  (k)  } 

1  ^  f.  \  p[zf^k)\>^;:f/,Ai-'{k),D^-'-)x' 

njii  P{Aj.‘:H7|A)-'(A;),D'=->} 

Then  we  can  easily  obtain  the  likelihood 

Aj‘{k)  &  p(7'(fc)iA77,Ai-‘W,B‘“‘) 

=  |2<rSf  (A)j'*  exp  |-1  (x”'  (k)  -  4{I=))'S]'‘  (k)-'  (xf  (k)  -  4  {k))  | 

where 

57^'  (k)  =  H  (k)  Pl^  (k)  H  (k)'  -  H  (k)  {k)  -  (fe)  H  (k)'  +  Rp  (fc) 

with 

Pt{k)  4  E[xpk)x1{ky} 

Rp{k)  4  E  [zp  {k)  zp  {k)'} 

Uyp{k)  4  E[xpk)zp{ky} 

Uppk)  4  E[f;^^{k)xpk)'] 

xf  {k)  =  xi  {k)  -  xf  (k) 
zp{k)  =  Hxi{k)-zp{k) 

and  x^  (k)  is  the  most  recent  fused  estimate/prediction  for  target  U  for  the  hypothesis 
being  considered. 


Hence  combining  the  cases  of  j  =  1,  1  <  j  <  J  and  again  assuming  that  no  two 
resolved  ground  tracks  which  are  due  to  the  same  target  can  be  associated  with  the  same 
propagation  path  gives  the  following  probability  update  for  the  hypothesis 
using  the  measurement  of  update  k,  5j  (k): 

p{>^,':~/\Ai(k),Dk-^}  = 

_ AT‘(k)P{K,\:r/  |A>--(A;),C*-1} _ 

1  Ml  Mj  _  _  ^  ' 

E  ■■  E  E  E  A”y*)p{AS,:7|Ai-i(*).£>‘-i} 

fSi=l  iij=l  mi=l  mj=l 
mj^Sj 

where,  for  j  =  1: 


=  Ar(/c)Ap(xrW|A7«7.-D''^‘) 

AfW  =  Af'{k)^p[zf‘{k)\>4;:f/,D>=-') 
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and  for  1  <  j  <  J: 

AfW  =  p(2rwiA?.'.:„7,A>-' (*),/?<=-■) 
kf(k)  A  p[zj•{k)\\t;::Z\^‘-'{k),D^-') 

The  likelihood  for  both  cases  is  then 

{k)  =  |27r5f^'  (/c)p^  exp  (fc)'  Sf^ 

where 

and 

^  E[x1{k)x^{kr] 

Rj^{k)  ^  E  {k)  {k)'] 

U’lp{k)  4  E{x'l{k)z';^{ky} 

xf  (A:)  =  Xi  (A:)  -  xf  (A:) 
z^^{k)  =  Hxi[k)  -  z'^^  {k) 

and  xf  (Ac)  is  the  most  recent  fused  estimate/prediction  for  target  ti  that  z^^  (Ac)  is  assumed 
to  be  associated  with  for  the  hypothesis  being  considered. 

Consider  now  the  estimate  x^"^^  (Ac)  of  Xj  (Ac)  and  its  covariance  (Ac)  obtained  by 
the  fusion  of  xk{k),h  =  —  1  with  z^^  (Ac).  This  is  again  derived  using  the 

fundamental  equations  of  linear  estimation  and  essentially  the  same  reasoning  as  shown 
earlier  for  Ac  =  0 .  The  resulting  fusion  equations  are 

(k)  =  x}  (k)  +  [pt  (k)  H  (*)'  -  {k})  S]‘>  (*)-‘  (2“'  (*)  -  2?  (*:))  (32) 

Pt*'  (k)  =  Pt  {k)  -  [Pt  {k)  H  (k)'  -  U^P‘  (k))  Sp  (k)-‘  (h  {k)P^  (k)  -  U^‘'‘  (k)) 

(33) 
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9  Fusion  of  Microwave  Radar  and  Multiple 

OTHR  tracks 

During  the  development  of  the  DMPTF  algorithm,  it  was  always  kept  in  mind  that 
the  algorithm,  if  possible,  should  be  extendable  to  multi-sensor  applications  thus  giving  a 
single  coherent  algorithm  for  fusing  multipath  OTHR  tracks  and  tracks  from  other  somrces 
such  as  microwave  radar  and  GPS.  This  aim  was  in  fact  achieved,  so  that  microwave  radar 
tracks  (or  GPS  tracks)  can  be  easily  incorporated  into  the  fusion  algorithm  by  treating 
estimates  from  the  microwave  radar  in  a  similar  manner  to  the  OTHR  tracks. 

Some  properties  of  the  microwave  radar  tracks  that  should  be  noted,  are 

•  The  tracks  generally  only  have  a  single  propagation  path  associated  with  them  (ex¬ 
cept  in  occasional  anomalous  circumstances),  hence  any  two  tracks  from  the  same 
microwave  radar  will  not,  in  general,  be  from  the  same  target. 

•  The  microwave  radars  are  generally  not  co-located  with  the  OTHR  hence  transfor¬ 
mations  need  to  be  provided  to  a  common  ground  coordinate  system. 

•  Covariance  matrices  for  the  microwave  track  estimates  are  often  not  available  in 
legacy  systems,  hence  approximations  of  these  covariance  matrices  may  need  to  be 
made. 

Let  us  now  consider  how  microwave  radar  tracks  would  be  incorporated  into  the  hy¬ 
pothesis  tree.  Essentially  the  microwave  radar  tracks  can  be  treated  almost  identically  to 
the  OTHR  tracks.  The  microwave  radar  tracks  arrive  asynchronously  with  respect  to  the 
OTHR  tracks;  however,  this  is  simply  accommodated  through  the  use  of  the  predictions 
that  are  produced  by  the  MPTF  target  model.  For  example,  if,  say,  two  OTHR  tracks 
result  in  the  creation  of  the  hypothesis  tree  at  update  k  —  0  and  a  microwave  radar  track 
is  introduced  at  time  k  =  1  then  the  tree  at  A:  ===  0  is  simply  extrapolated  to  A;  =  1  and  then 
extended  to  3  tracks  using  the  microwave  radar  data.  Figure  10  shows  the  path  indepen¬ 
dent  hypothesis  tree  for  this  example.  The  path  dependent  tree  will  of  course  incorporate 
the  number  of  paths  occuring  with  each  of  the  OTHR  tracks  and  only  one  “path”  for  the 
microwave  radar  track.  The  “path”  for  the  microwave  radar  track  is  the  transformation 
from  the  microwave  radar’s  coordinates  to  the  common  ground  coordinates  that  are  used. 

The  code  for  incorporating  microwave  radar  tracks  into  the  DMPTF  algorithm  has 
been  implemented  in  the  test-bed,  and  is  expected  to  be  tested  in  the  near  future  on  real 
track  data.  The  main  extension  that  was  required  was  to  derive  the  transformations  for  the 
microwave  radar  tracks  to  the  common  ground  coordinates.  The  JFAS  ground  coordinates 
were  chosen  to  be  the  common  coordinate  system,  ie,  range  and  azimuth  relative  to  the 
receive  array,  with  the  receive  array  boresight  corresponding  to  zero  azimuth. 

The  same  principles  that  have  been  applied  to  fusing  the  OTHR  multipath  tracks  with 
microwave  radar  tracks  apply  to  fusion  of  OTHR  and  GPS  tracks,  as  well  as  fusion  of 
multipath  tracks  from  multiple  OTHRs  in  an  overlapping  network.  With  GPS  tracks  the 
method  of  fusion  with  OTHR  tracks  is  essentially  identical  to  that  for  fusion  of  microwave 
radar  and  OTHR  tracks.  In  the  case  of  multiple  OTHRs,  each  OTHR  will  have  a  different 
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Figure  10:  Path  independent  hypothesis  tree  for  two  OTHR  tracks  at  time  k  =  0  and  one 
microwave  radar  track  at  time  k  —  1 

set  of  multiple  paths,  but  this  is  easily  accommodated  in  the  DMPTF  algorithm.  With 
regard  to  coordinate  systems,  a  common  ground  coordinate  system  will  be  required  as 
opposed  to  the  polar  coordinates  referenced  to  JFAS  that  is  used  at  present.  In  principle 
this  poses  no  difficulties  other  than  the  requirement  for  more  computation  to  be  performed 
to  convert  to  the  common  coordinates. 

A  point  worth  making  here  is  that,  in  general,  microwave  radar  tracks  and  GPS  tracks 
are  much  more  accurate  than  OTHR  tracks.  This  fact  can  be  used  to  improve  coordinate 
registration  (CR)  of  the  OTHR  where  microwave  radar  or  GPS  tracks  are  available.  For 
example,  if  the  registration  of  the  OTHR  is  initially  good  enough  to  perform  correct 
association  of  a  group  of  multipath  OTHR  tracks  with,  say,  their  corresponding  microwave 
radar  track,  then  the  offset  between  the  OTHR  tracks  and  the  microwave  radar  track  can 
be  used  to  perform  corrections  to  the  OTHRs  slant  to  ground  coordinate  transformations. 
Because  the  height  of  ionospheric  layers  can  be  expected  to  not  change  abruptly,  the 
corrections  will  be  able  to  be  utilized  for  regions  around  the  tracks  as  well.  The  extent  of 
these  regions,  and  the  degradation  in  the  accuracy  of  the  CR  corrections  as  a  function  of 
distance  from  the  registered  tracks  would  need  to  be  determined  by  experiment. 
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10  Summary  and  Conclusions 

A  dynamic  multipath  fusion  (DMPTF)  algorithm  which  is  capable  of  fusing  over-the- 
horizon  radar  (OTHR)  multipath  tracks  and  non-OTHR  tracks  (e.g.  microwave  radar  or 
GPS),  as  well  as  dealing  with  multipath  tracks  from  OTHR  networks  is  presented  in  this 
report.  The  algorithm  achieves  this  through  a  very  general  model  based  approach  which 
can  deal  with  both  the  multipath  effects  within  a  single  OTHR  as  well  as  asynchronicity 
between  multiple  sensors.  Important  advances  of  the  algorithm  presented  in  this  report 
over  earlier  related  work,  ie,  [17],  [18]  and  [19]  are: 

•  The  new  algorithm  can  associate  and  fuse  asynchronous  track  data,  thus  making 
combined  multisensor-multipath  track  fusion  possible;  particularly  of  interest  is  the 
algorithm’s  direct  applicability  to  fusion  with  microwave  radar,  GPS  and  multiple 
OTHR  tracks. 

•  The  new  algorithm  takes  into  account  temporal  relationships  between  multipath 
and  multisensor  tracks  as  well  as  the  instantaneous  spatial  relationships.  This  is  of 
theoretical  importance  in  addition  to  benefitting  performance. 

•  An  adjunct  algorithm  which  can  effectively  deal  with  increases  and  decreases  in  the 
number  of  tracks  with  time  in  DMPTF’s  hypothesis  tree  has  been  developed.  This 
is  important  for  achieving  temporally  consistent  association  hypotheses  and  fused 
tracks.  Without  this,  each  time  a  multipath  track  drops  out,  a  new  hypothesis  tree 
would  need  to  be  constructed,  losing  all  the  information  in  the  previous  tree. 

•  An  effective  and  efficient  pruning  algorithm,  as  well  as  a  clustering  algorithm  have 
been  developed  which  together  enable  association  and  fusion  of  very  large  numbers 
of  tracks.  This  has  been  tested  on  the  JFAS  OTHR  in  an  operational  setting  showing 
that  the  DMPTF  algorithm  together  with  the  associated  pruning  and  clustering  can 
successfully  deal  with  all  the  tracks  in  the  entire  coverage  area  of  the  radar  in  real 
time. 

•  A  set  of  equations  have  been  derived  for  the  association  and  fusion  of  the  measure¬ 
ments  from  which  the  multipath  and  multisensor  tracks  are  created.  This  enables 
an  alternative  fusion  approach  where  some  sources  of  track  dependence  are  avoided. 

The  DMPTF  algorithm  has  been  implemented  in  C-h-f  in  a  test-bed,  for  trialling  with 
simulated  and  real  data,  and  also  implemented  as  a  prototype  on  the  JFAS  OTHR  for 
testing  in  an  operational  setting  and  for  collection  of  data  for  statistical  analysis.  To  date, 
testing  has  centred  on  multipath  data  from  the  JFAS  radar.  Code  has  also  been  imple¬ 
mented  for  fusion  with  microwave  radar  tracks  which  awaits  testing.  Considerable  eflFort 
was  invested  in  the  development  of  the  C-f— (-  code  for  DMPTF  to  achieve  modularity,  ex¬ 
pandability,  and  maintainability,  and  a  semi-automated  software  documentation  program 
was  implemented.  As  an  outcome,  it  should  be  relatively  easy  to  further  develop  the  code 
and  transfer  it  to  another  target  system  such  as  JORN  if  required. 

With  regard  to  performance,  testing  of  DMPTF  with  OTHR  multipath  tracks  (ie,  real 
data)  has  been  performed  on  a  development  test-bed  from  early  stages  in  its  development. 
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Recently  an  extended  trial  of  the  algorithm  has  been  performed  on  the  JFAS  radar,  and 
a  large  amount  of  data  collected  for  statistical  analysis.  While  the  statistical  analysis 
has  not  yet  been  completed,  early  indications  are  that  the  algorithm  performs  well  when 
good  coordinate  registration  (CR)  data  is  available.  However,  when  the  CR  data  is  poor, 
incorrect  associations  often  occur.  Hence,  at  the  present  level  of  development  of  DMPTF 
and  the  CR  system,  the  DMPTF  algorithm  is  capable  of  being  used  in  a  semi-automated 
fashion  to  give  advice  to  an  operator,  but  operator  overseeing  is  still  required.  Higher 
levels  of  automation  can  be  expected  with  further  development  of  both  DMPTF  and  CR. 

Further  operational  assessment  and  development  is  strongly  recommended.  In  partic- 

The  code  for  automatically  associating  and  fusing  multipath  OTHR  tracks  should 
be  further  operationally  assessed  and  enhanced. 

The  code  for  fusion  of  OTHR  and  microwave  radar  tracks  which  has  been  imple¬ 
mented  should  be  assessed. 

The  fusion  of  tracks  in  OTHR  networks  in  regions  of  overlapping  surveillance  cover¬ 
age  should  be  implemented  and  assessed. 

The  potential  for  improving  OTHR  coordinate  registration  accuracy  via  feedback 
from  other  surveillance  sensors  should  be  developed. 

Further  research  on  track  dependence  and  development  of  a  realistic  ionospheric 
model  (the  two  are  related)  should  bear  fruit  with  regard  to  improving  performance. 


ular: 


54 


DSTO-RR-0223 


Authors’  Contributions 

The  contributions  of  the  authors  were  as  follows: 

Peter  Sarunic  performed  the  first  principles  derivation  of  the  static  MPTF,  and  pro¬ 
posed  and  led  the  development  of  the  new  pruning  algorithm.  He  conceived  and  performed 
the  mathematical  derivations  of  the  dynamic  multipath  track  fusion  algorithm,  proposed 
the  extension  of  DMPTF  to  perform  fusion  of  selected  measurements  to  reduce  data  de¬ 
pendence,  and  derived  the  equations  for  doing  this.  He  also  wrote  the  core  of  the  C-|— h 
code  for  the  fusion  algorithms  and  led  the  subsequent  coding  and  testing. 

Kruger  White,  who  had  been  working  on  multipath  fusion  for  some  time  prior  to 
the  work  presented  in  this  report,  provided  much  useful  information  and  advice  based 
on  his  earlier  work  on  the  problem,  including  advice  on  fimitations  he  had  encountered 
with  the  earlier  fusion  algorithm.  He  also  developed  coordinate  transformation  software, 
participated  in  literature  surveys,  particularly  on  track  dependence,  and  proposed  the 
investigation  of  the  Covariance  Intersection  technique. 

Mark  Rutten  joined  the  team  shortly  after  the  commencement  of  the  work  which  is 
the  subject  of  this  report.  After  Peter’s  original  coding  of  the  static  MPTF  code  in  C-l— b, 
Mark  coded  extensions  to  include  pruning,  target  dynamics  and  clustering.  He  contributed 
to  the  development  of  the  K-best  search  algorithm  and  developed  the  theoretical  aspects 
of  clustering  (not  covered  in  this  report).  He  also  contributed  to  testing  of  algorithms. 


Acknowledgements 

The  authors  would  like  to  thank  Peter  den  Hartog  and  Kim  Kieu  for  their  efforts  in  the 
implementation  and  testing  of  the  DMPTF  software.  Peter  den  Hartog  performed  integra¬ 
tion  of  the  software  with  the  JFAS  radar,  organised  trials,  and  did  some  of  the  graphical 
display  work  for  the  test-bed  software.  He  also  developed  a  web-site  for  presentation  of 
algorithm  developments.  Kim  did  the  major  part  of  the  display  graphics  for  the  test-bed 
and  the  integration  of  DMPTF  with  the  graphics  software.  Both  Kim  Kieu  and  Peter  den 
Hartog,  contributed  to  the  testing  of  the  fusion  algorithms  and  associated  software.  Their 
efforts  are  much  appreciated.  While  John  Percival  did  not  participate  significantly  in  the 
work  presented  in  this  report,  he  played  a  key  role  in  the  earlier  work  on  Multipath  Tback 
Fusion  which  is  referenced  in  this  report,  for  which  his  efforts  are  recognized. 


55 


DSTO-RR-0223 


References 

1.  G.  Adelson-Velskii  and  E.  Landis.  “On  an  Information  Organization  Algorithm”,  Dok- 
lady  Akademica  Nauk  SSSR,  Vol.  146,  pp.  263-266,  1962. 

2.  M.  Azmoodeh,  “Abstract  Data  Types  and  Algorithms”,  MacMillan,  1988. 

3.  Y.  Bar-Shalom  and  T.E.  Fortmann,  “Tracking  and  Data  Association”,  Academic 
Press,  Inc.,  1988. 

4.  Yaakov  Bar-Shalom  and  Xiao-Rong  Li,  “Estimation  and  Tracking;  Principles,  Tech¬ 
niques  and  Software”,  Artech  House,  1993. 

5.  Yaakov  Bar-Shalom  and  Xiao-Rong  Li,  “Multitarget-Multisensor  Tracking:  Principles 
and  Techniques”,  Yaakov  Bar-Shalom,  1995. 

6.  L.  Comtet,  “Advanced  Combinatorics  -The  Art  of  Finite  and  Infinite  Expansions”,  D. 
Reidel  Publishing  Company,  Dordrecht,  Holland,  1974. 

7.  C.Y.  Chong,  S.  Mori,  and  K.C.  Chang,  “Distributed  multitarget  multisensor  track¬ 
ing”  in  Multitarget  Multisensor  Tracking:  Advanced  applications,  Y.  Bar-Shalom,  Ed., 
Norwood  MA:  Artech  House,  1990,  pp  247-295. 

8.  I.W.  Dali  and  D.J.  Kewley,  “Track  Association  in  the  Presence  of  Multi-Mode  Prop¬ 
agation”,  Proc.  International  Conference  on  Radar  92,  IEEE  Publication  365,  pp. 
70-73,  1992. 

9.  C.A.R.  Hoare,  “Proof  of  a  Program:  FIND”,  CACM,  Vol  14,  No.l,  pp  39-45,  Jan 
1971. 

10.  J.H.  Kingston,  “Algorithms  and  Data  Structures”,  Addison  Wesley,  1990. 

11.  J.L.  Krolik  and  R.H.  Anderson,  “Maximum  Likelihood  Coordinate  Registration  for 
Over-the-Horizon  Radar”,  IEEE  Trans,  on  Signal  Proc.  45,  935-959,  1997. 

12.  T.  Kurien,  D.  Logan  and  W.P.  Berry,  “Fusion  of  OTH  Radar  Data”,  Proc.  Seventh 
Joint  Service  Data  Fusion  Symposium,  pp.  221-235,  1994. 

13.  M.E.  Liggins,  C.Y.  Chong,  1.  Kadar,  M.G.  Alford,  V.  Vannicola  and  S.  Thomopoulos, 
“Distributed  Fusion  Architectures  and  Algorithms  for  Target  Tracking” ,  Proc.  of  the 
IEEE  85,  1,  pp  95-107,  1997. 

14.  S.  Mori,  K.A.  Demetri,  W.H.  Barker  and  R.N.  Lineback,  “A  Theoretical  Foundation 
of  Data  Fusion  -  Generic  Track  Association  Metric” ,  Proc.  Seventh  Joint  Service  Data 
Fusion  Symposium,  pp.  585-594,  1994. 

15.  Lucy  Y.  Pao,  “Measurement  Reconstruction  Approach  for  Distributed  Multisensor 
Fusion”,  Journal  of  Guidance,  Control,  and  Dynamics,  Vol  19,  No.  4,  pp.  842-847, 
Jul-Aug  1996. 

16.  Lucy  Y.  Pao  and  Michael  Kalandros,  “Algorithms  for  a  Class  of  Distributed  Architec¬ 
ture  Tracking” ,  Proc.  of  the  American  Control  Conference,  Albuquerque,  New  Mexico, 
pp.  1434-  1438,  June  1997. 


56 


DSTO-RR-0223 


17.  D.J.  Percival  and  K.A.B.  White,  “Multipath  Track  Fusion  for  Over-the-Horizon 
Radar”,  Proc.  SPIE  3163,  pp.  363-374,  1997. 

18.  D.J.  Percival  and  K.A.B.  White,  “Multipath  Coordinate  Registration  and  TVack 
Fusion  for  Over-the-Horizon  Radar”,  accepted.  Signal  Proc.  Special  Issue  for  the 
DSTO/ AFOSR  Signal  Processing  Workshop,  Victor  Harbour,  South  Australia,  1997. 

19.  D.J.  Percival  and  K.A.B.  White,  “Multihypothesis  Fusion  of  Multipath  Over-the- 
Horizon  Radar  TVacks”,  Proc.  SPIE  3373,  pp.440-451,  1998. 

20.  R.L.  Popp,  K.R.  Pattipati  and  Y.  Bar-Shalom,  “m-best  SD  Assignment  Algorithm 
with  Application  to  Multitarget  Tracking” ,  Proc.  SPIE  3373,  pp.475-495,  1998. 

21.  W.  Pugh,  “Skip  Lists:  A  Probabilistic  Alternative  to  Balanced  Trees,”  Communica¬ 
tions  of  the  ACM,  Vol.  33,  pp.  668-676,  June  1990. 

22.  G.W.  Pulford  and  R.J.  Evans,  “A  Multipath  Data  Association  Ttacker  for  Over- 
the-Horizon  Radar”,  IEEE  Trans,  on  Aerospace  and  Electronic  Systems,  Vol.  34,  pp 
1165-1183,  October  1998. 

23.  D.B.  Reid,  “An  Algorithm  for  Tracking  Multiple  Targets”,  IEEE  Trans.  on  Automatic 
Control,  December  1979. 

24.  P.W.  Sarunic,  M.  Rutten,  D.J.  Percival  and  K.A.B.  White,  “Hypothesis  Management 
for  Over-the-Horizon  Radar  Multipath  Track  Fusion”,  Proc.  EuroFusion99  Interna¬ 
tional  Conference  on  Data  Fusion,  pp  217-224,  Oct.  1999. 

25.  P.W.  Sarunic,  M.G.  Rutten,  ’’Over-the-Horizon  Radar  Multipath  Track  Fusion  Incor¬ 
porating  Track  History” ,  Proc.  Pusion2000  Conference,  pp  TuCl-13  to  TuCl-19,  July 
2000. 

26.  R.  Sedgewick,  “Algorithms” ,  Second  Edition,  Addison-  Wesley,  1988. 

27.  J.K.  Uhlmann,  “General  Data  Fusion  for  Estimates  with  Unknown  Cross  Covari¬ 
ances”,  Proc.  SPIE  2755,  pp.  536-547,  1996. 

28.  J.  Zhu,  R.E.  Bogner,  A.  Bouzerdoum  and  M.L.  Southcott,  “Application  of  Neural 
Networks  to  Track  Association  in  Over  the  Horizon  Radar”,  Proc.  SPIE  2233,  pp. 
224-235,  1994. 


57 


DSTO-RR-0223 


Appendix  A:  Summary  of  Errors  in  Previously 
Derived  MPTF  Equations 

This  appendix  summarizes  errors  found  in  some  of  the  static  equations  presented  in  [17] 
and  [19].  The  equations  were  determined  by  analogy  arguments  from  multisensor  fusion 
work  published  in  the  literature;  however,  the  translation  to  the  multi-path  fusion  prob¬ 
lem  domain  incurred  some  mistakes.  The  equations  in  question  are  the  counterparts  of 
equations  2,  3,  and  5  in  section  2  of  this  report.  In  [17]  and  [19]  two  probability  update 
equations  are  presented.  The  equations  are  presented  below  using  the  notation  of  [19]. 
The  reader  is  referred  to  the  papers  in  question  for  a  detailed  description  of  the  notation, 
as  well  as  a  description  of  how  the  equations  were  obtained. 

Pr{AlZ''UZ^,D''}  =  (Al) 

I  C-lPr{A|Z^T>''}  e  A,  ^  A, 

\  C-lPr{A|Z^iA^}  ,D^),  =  f U  G  A 

The  first  equation  is  the  probability  update  for  a  track  estimate  which  is  assumed  to  be 
due  to  a  new  target  for  the  hypothesis  in  question.  The  second  is  the  probability  update 
for  a  track  estimate  that  represents  the  same  target  as  at  least  one  previously  hypothesized 
target. 

The  second  of  equations  Al  has  a  likelihood  term  which  in  reference  [17]  is  obtained 
using  the  following  relationship 

^  U  Zf  U  Z^)  =  (A2) 

27r  -f  {Pir^  +  PiT^'^) 

whereas  in  [19]  the  likelihood  term  is  determined  using  the  equation  below 

^^-1  |27rn|-2exp  (A3) 

with 

n  = 

and  where  p  is  the  a  priori  target  density. 

Firstly,  in  both  of  the  probability  update  equations  Al  there  is  a  normalization  constant 
C"^  The  way  to  determine  its  value  is  not  given  in  either  of  the  papers;  however,  in  an 
implementation  of  the  static  MPTF  algorithm  based  on  the  above  papers,  C  was  taken  to 
be  the  sum  of  the  unnormalized  probabilities  at  level  k  of  the  hypothesis  tree,  ie, 

c=  Pr„ma/i.ed{A|Z''UZ^P'^}  (A4) 

All  A  at 
level  k 
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This  is,  in  fact,  incorrect  as  can  be  seen  from  equation  2  in  section  2  of  this  report. 

Secondly,  comparing  the  first  equation  with  equations  2,  3  in  section  2,  shows  that  it 
is  missing  a  1/1^  term.  This  is  of  course  important,  as  it  affects  the  relative  magnitude 
of  the  probabilities  for  new-target  hypotheses  when  compared  to  hypotheses  assuming  an 
already  existing  target. 

Thirdly,  because  equation  A4  was  mistakenly  used  for  normalization,  it  made  the 
second  probability  equation  of  equations  Al  dimensionally  inconsistent  if  the  likelihood 
equation  of  reference  [17]  (ie,  equation  A2)  is  used  for  the  likelihood  term.  To  compensate 
for  this,  an  adjustment  was  made  to  the  likelihood  equation  (ie,  adding  the  p~^  term, 
where  p  is  an  assumed  target  density),  giving  equation  A3  which  was  presented  in  [19]. 
Unfortunately,  this  adjustment  still  did  not  correct  the  equations,  as  can  be  seen  by 
comparison  with  equations  2  and  5  in  section  2  of  this  report. 
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